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Abstract 

Achievements in the heavy quark theory over the last decade are reviewed, with 
the main emphasis put on dynamical methods which quantify nonperturbative ef- 
fects via application of the Operator Product Expansion. These include the total 
weak decay rates of heavy flavor hadrons and nonperturbative corrections to heavy 
quark sum rules. Two new exact super convergent sum rules are derived; they differ 
from the known ones in that they are finite and normalization point independent 
in perturbation theory. A new hadronic parameter E is introduced which is a spin- 
nonsinglet analogue of A = Mb — m&; it is expected to be about 0.25 GeV. The first 
sum rule implies the bound p 2 > 3/4 for the slope of the Isgur-Wise function. The 
heavy quark potential is discussed and its connection to the infrared contributions 
in the heavy quark mass. Among applications extraction of \V c b\ from the total 
semileptonic and from the B^D* zero recoil rates is addressed, as well as extract- 
ing \V u b\ from r s i(6— >u). Practical aspects of local quark-hadron duality are briefly 
discussed. 
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1 Introduction 



Quark-gluon dynamics are governed by the QCD Lagrangian 

C = - 1 -G%G% + J2^P-m^ (1) 

where G a is the gluon field strength tensor, and ip are quark fields. Masses of u and 
d quarks are only a few MeV, much smaller than Aqcd, and in most applications 
of hadron physics u and d can be considered as massless. The strange quark mass 
is about 150 MeV, only a little smaller than Aqcd- Nevertheless, there is ample 
evidence that treating s quark as light is justified, and corrections to the (light) 
SU(3) symmetry are reasonably small. Thus, m s lies essentially below the actual 
typical hadronic QCD scale /ihadr ~ 500 to 700 MeV. 

Those quarks Q for which uiq ^> /ihadr ~ (2 -3) Aqcd are heavy quarks. The 
sixth t quark is the heaviest, m t ~ 170 GeV. However, it is too heavy. Its width 
due to the 'semiweak' decay t — > b + W + is T t ss 1 GeV. It decays too fast for the 
t-hadrons, the bound states with light quarks to be formed. 

The best candidates for application of the Heavy Quark Expansion (HQE) are 
beauty hadrons. b quark is heavy enough to confidently use the expansion in l/m b , 
yet the leading corrections to the heavy quark limit are not negligible. 

The charmed quark c can be called heavy only with some reservations. While 
in some cases this yields a reasonable approximation, m c often appears manifestly 
too low for a quantitative treatment of charmed physics in the 1/toq expansion. 
There is no universal answer here, and considerable care must be exercised in every 
particular case. 

Usually considered heavy flavor hadrons are composed from one heavy quark 
Q, and a light cloud: a light antiquark q, or diquark qq, together with the gluon 
medium and light quark- antiquark pairs. The role of the gluon medium is to keep all 
'valence' constituents together, in a colorless bound state which will be generically 
denoted by Hq. Therefore, we have the following simplified picture of a heavy flavor 
hadron. The heavy quark has a small size ~ 1/m.Q, and is surrounded by a static 
Coulomb-like color field Aq at small distances. Non-Abelian selfinteraction slightly 
modifies the potential, but the nonlinearity is driven by the coupling a s /V and 
is not significant. At larger distances the selfinteraction strengthens, at R > Aq CD it 
is completely nonperturbative: the soft modes of the light fields are strongly coupled 
and strongly fluctuate. 

In weak decays the standard situation is that the external (to QCD) forces like 
W bosons interact with the heavy quark, say, instantaneously replace the b quark 
by c quark generally changing its velocity. Such an event excites first the typical 
modes of the heavy quark, hard gluons with k typ ~ toq. Since 

a s (hyp) ~ a s (m Q ) < 1 
perturbation theory is adequate there. 
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An actual decay process, nevertheless, eventually runs into the strong-interaction 
nonperturbative domain of k ~ uo ~ Aqcdj where u denotes the characteristic fre- 
quencies. The final hadronization dynamics shaping the hadrons observed in experi- 
ment, is a result of soft nonperturbative physics which is responsible for confinement. 
Since the complicated final state dynamics involve u <C rriQ, to deal with nonper- 
turbative effects one can use the nonrelativistic expansion for the heavy quark. 

The main subject of the HQE is nonperturbative physics. First, the perturbative 
corrections are conceptually simple, even if the actual computations are often cum- 
bersome and, beyond the first-order effects, require sophisticated state-of-the-art 
technique. 

Second, the perturbative corrections are calculated in the full QCD rather than 
in the effective low-energy theory, since they come just from the gluon momenta k ~ 
iriQ where the nonrelativistic approximation is not applicable. Still, the interplay of 
the perturbative and nonperturbative effects is quite nontrivial and often involves 
theoretical subtleties. 

The treatment of the nonperturbative effects is a nontrivial problem, and differ- 
ent methods of QCD are used here. The basic tool for all of them in heavy quarks 
is the Wilson operator product expansion (OPE) |J. 

1.1 Nonrelativistic expansion 

The main simplification of a nonrelativistic treatment is that the number of heavy 
quarks nq and antiquarks ftq are separately conserved. Propagation with usual 
relativistic Green function for the heavy quark contains also the process of the 
QQ pair creation if the time ordering of the vertices along the line is reversed 
somewhere. In the nonrelativistic kinematics, however, such configurations yield 
a power- suppressed contribution, since the virtuality of the intermediate state (the 
energy denominator, in the language of noncovariant perturbation theory of quan- 
tum mechanics) is of the order of 2mg. In order to observe such processes as real 
ones, it would be necessary to supply energy at least as large as u ~ 2/bq. 

Heavy quarks can appear also in closed loops, for example, in the processes of 
the virtual gluon conversion. For heavy quarks such effects are also suppressed, 
~ k 2 jm 2 Q if the gluon momentum k is much smaller than mq. 

As a result, the field-theoretic, or second-quantized description of the heavy 
quark becomes redundant, and it is sufficient to resort to its usual quantum-mecha- 
nical (QM) treatment. For example, the wavefunction of a heavy flavor hadron takes 
the form 

*a [% {flight}] , (2) 

where xq is the heavy quark coordinate, xu g ht generically represents an infinite 
number of light degrees of freedom; index a describes the heavy quark spin. 

A relativistic S — ^ particle has four components, i.e. a = 1...4. The non- 
relativistic spinor ty(x) has only two of them, a = 1,2 describing two spin states. 
The remaining a = {3, 4} components describe antiparticles which decouple in the 
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nonrelativistic theory: 



Q( x ) = I J J ) x ( x )^^^n ( 3 ) 




The nonrelativistic Hamiltonian of the spin-| particle is the well known Pauli Hamil- 
tonian:Fl 

Hpauli = -A) + V - ^ + ^ • (4) 

The last operator is written with the coefficient appearing for an 'elementary' point- 
like particle. Presence of the interaction generally renormalizes its strength, the 
chromomagnetic moment of the heavy quark. In the heavy quark limit mg — > oo 
only the first term survives while the last two terms describing space propagation 
and interaction of spin with the chromomagnetic field B, disappear. An infinitely 
heavy quark is static and interacts only with the color Coulomb potential.0 In turn, 
the heavy quark is the source of a static color field independent of the heavy quark 
spin. The actual dynamics of the heavy quark spin reveal itself in full only when 
the quark interacts with quanta having large momentum \q\ > w,q. 

It is easy to illustrate the b quark spin-independence of the strong forces in the 
following way. One can imagine the QCD world where, instead of the actual b quark 
there exists a scalar spinless b quark with the same mass. The usual spin-0 B meson, 
bq would become a spin-| particle B ~ bq. The A& baryon bud having spin-| would, 
in turn become a scalar spinless bud state; would have spin 1. Nevertheless, at 
rrib — > oo the properties of B and B or A& and A& would be identical. The two 
degenerate spin states of B are counterparts of B and B* of the actual QCD. 

An immediate question comes to one's mind at this point: what about the 
relation between spin and statistics? In such a gedanken operation we replace spin- 
\ hadrons by S = ones, and vice versa, i.e. interchange fermions and bosons. 
It is clear, however, that for the states or processes with a single heavy quark the 
statistics symmetry properties do not play a role. 

Such an independence of the strong dynamics of the heavy quark spin is called 
the heavy quark Spin Symmetry. 

The rest mass itiq of the heavy quark also enters in a trivial way: the Hamiltonian 
simply contains (nQ+nQ)rriQ. Since both are fixed, it is an overall additive constant. 
For a moving quark this constant is E = <J ttiq + p 2 = mg^l + v 2 . Therefore, the 
actual dynamics is not affected by the concrete value of the mass mq. This is the 
Heavy Flavor Symmetry which states, for example, the equal properties of charmed 
and beauty hadrons to the extent they both can be considered heavy enough. 

The heavy flavor symmetry leads also to certain scaling behavior of the transition 
amplitudes with heavy flavor hadrons: the amplitudes depend on their velocities 

1 We use the convention where the coupling g s is absorbed in the gauge fields; this is convenient 
for nonperturbative analysis. 

2 Speaking of the Coulomb interaction in QCD we mean the interaction with the timelike com- 
ponent Aq of the color gauge potential. It differs from the simple 1/R electrodynamic potential. 
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rather than on the absolute values of momenta: 

(pin pout \ 

- 9 -,pS b ; — ,ph • (5) 
m Q m Q > J 

Here P denote the momenta of the heavy flavor hadrons and p refer to other partic- 
ipating particles. Such a scaling behavior is valid only with respect to the soft part 
of the interaction, and no hard gluons with k ~ mq are involved. 

The one-particle (or QM) description of the heavy quark degrees of freedom is 
the key simplification of the nonrelativistic expansion. The light cloud, however, 
still requires a full-fledged field-theoretic treatment. Even considering the static 
limit niQ — ► oo where only interaction with the Coulomb field A$ remains, the 
latter strongly fluctuates, in contrast to simple potential QM models where the 
corresponding potential V(x) is a c-number function of coordinates. 

For this reason even the quantum mechanics of heavy quarks is highly nontriv- 
ial. Exploiting the symmetry properties of the heavy quark interactions does not 
require understanding of these complicated strong interaction dynamics. This was 
the main field of applications at an early stage of theoretical development of heavy 
quark physics. The recent progress is mainly related to a better treatment of basic 
properties of this complicated strongly interacting system via application of dynamic 
QCD methods based on the short distance expansion. 

1.2 Operator Product Expansion 

The basic theoretical tool of the heavy quark theory is the Wilson operator product 
expansion 0. The idea of the OPE was formulated by K. Wilson in the late 60's, 
originally in the context of the statistical problems which are closely related to 
the field theories in Euclidean space. This idea, in general, is a separation of effects 
originating at large and small distances. The application to real physical processes in 
Minkowski space is often less transparent and technically more complicated, however 
it is always based on the same underlying concept. 
The original QCD Lagrangian 

c = -\G; u G%+J2q(^-™ q )q+EQ(^- m Q)Q = c light +J2Q(^-m Q )Q (6) 

9 Q Q 

is formulated at very short distances, or, equivalently, at a very high normalization 
point p = Mo, where M is the mass of an ultraviolet regulator. In other words, 
the normalization point is assumed to be much higher than all mass scales in the 
theory, in particular, p ^> mq. An effective theory for describing the low energy 
properties is obtained by evolving the Lagrangian from the high scale Mq down to a 
lower normalization point p. This means that we integrate out, step by step, all high 
frequency modes in the theory thus calculating the Lagrangian C(p). The latter is a 
full-fledged Lagrangian with respect to the soft modes with characteristic frequencies 
less than p. The hard (high frequency) modes determine the coefficient functions 
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in C(fx), while the contribution of the soft modes is hidden in the matrix elements 
of (an infinite set of) operators appearing in C(fx). The value of this approach, 
outlined by Wilson long ago [Q, has become widely recognized and exploited in 
countless applications. The peculiarity of the heavy quark theory lies in the fact 
that the in and out states contain heavy quarks. Although we integrate out the 
field fluctuations with the frequencies down to fx <C tjiq, the heavy quark fields 
themselves are not integrated out since we consider the sector with nonzero heavy 
flavor charge. The effective Lagrangian £(fx) acts in this sector. Since the heavy 
quarks are neither produced nor annihilated, any sector with the given Q-number is 
treated separately from all others. 

2 Basics of the Heavy Quark Theory 
2.1 Effective Hamiltonian 

The strategy for integrating out virtual degrees of freedom to pass on to an effective 
theory of heavy particles is described in the textbooks. The nonrelativistic fermion 
field Q(x) has four degrees of freedom. Two of them, in Eq. (§) are nearly 

on-shell and two xi x ) are highly virtual describing excitation of antiquarks. One 
needs to integrate out first the antiparticle fields x{ x )- F° r simplicity, we consider 
this in the rest frame. 

Integrating out the antiparticle fields can be done straightforwardly since the 
QCD Lagrangian is bilinear in the quark fields. In this way one would obtain the 
(tree level version of the) so-called Lagrangian of HQET £hqet- It is a correct 
nonrelativistic Lagrangian up to the first order in 1/mq. 

This does not complete the program, however: one yet has a full 'particle' field 
\I/(x) which includes all frequencies from to oo. Such a problem does not show up 
in the potential models where modes with k 3> Aqcd axe not excited. However, in 
actual QCD all radiative corrections would diverge due to large momentum gluons. 

Therefore, besides the antiparticle fields one needs to integrate out hard gluons 
and the high frequency components of the heavy quark field *f>(x) itself, those for 

— * 

which \k\,u > fx. The scale fx is the normalization point of the effective theory. 
Since the configurations we integrate out depend on fx, the remaining effective low 
energy theory is also /i-dependent. 

In practice, we want to have fx <C mj, and, actually, as low as possible. On the 
other hand, fx must still belong to the perturbative domain. In practice this means 
that the best choice routinely adopted for applications is fx ~ a few x Aqcd, from 
0.7 to 1 GeV. All coefficients in the effective Lagrangian obtained in this way are 
functions of fx. 

To summarize, in treating heavy quarks we separate all strong interaction effects 
into 'hard' and 'soft' introducing a normalization scale fx. To calculate the effect 
of the short distance (perturbative) physics we use the original QCD Lagrangian 
Eq. (||). The soft physics is treated by the nonrelativistic Lagrangian where the 
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heavy quarks are represented by the corresponding nonrelativistic fields (fiQ-. 
£eff = -j(^ + ^2q{ip-Tri q )q + ^2{-mQ(p^(pQ + (p^iD (pQ - 

9 Q 



oiD 
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2mQ ^ ' 8m, 

where 



(DE) + a-{Exn-nxE} (p Q \ (7) 



n = iD = p- A , (aifff = (an) 2 = n 2 + SB ; (8) 

for simplicity only the tree level //-independent coefficients are given. By the stan- 
dard rules one constructs from the Lagrangian (|7]) the corresponding Hamiltonian. 
The heavy quark part takes the form f\ 

H Q = -A + -^—{if 2 + aB) + — -z- \-(DE)+a-{Ex7f-7fxE}] +0(l/m%) . 
2m Q 8m Q L J 

(9) 

The first term in the l/m^ part is called the Darwin term and the second one is the 
convection current (spin-orbital, or LS) interaction. 

Since the external interactions (electromagnetic, weak etc.) are given in terms 
of the full QCD fields Q(x), one needs also the relation between Q(x) and ^Pq{x): 



Let us briefly recall the textbook procedure for obtaining the nonrelativistic 
Lagrangian. One starts with 

Cl avy = Q(x)(iP-m Q )Q(x) (11) 

and factors out of the Q(x) field the "mechanical" time- dependent factor associated 
with the rest energy mq\ 

Q(x) = e- im ^Q(x) . (12) 
In an arbitrary frame moving with four-velocity t> M it takes the following form: 

Q( x ) = e - im Q v ^Q{x) . (13) 

Then 

iDpQ(x) = e {m Q v, + tt m ) Q(x) , tt m = P M - m Q v, . (14) 

The Dirac equation (ip — rriQ) Q = takes the form (now the tilde on Q is omitted) 

^Q = ^Q, *«Q = -^p^Q. (15) 

2 2mr> 2mn 



3 The leading term Aq is often omitted here. It is then implied that the time evolution operator 
is 7To = i-^ + Aq rather than in the usual Schrodinger equation. This can be consistently carried 
out through the analysis. 
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which allows one to exclude the small low components —^-Q[x) expressing them via 
the 'large' upper components -^-Q(x). For example, we have the useful identity 

QQ = QloQ + 2QN2) Q = QloQ + Q \ Q + total derivative . (16) 
V > 2m z Q 

A subtlety emerges on this route that must be treated properly: at order 1/ttVq 
the time derivative do appears with the nontrivial coefficient depending, for example, 
on the gluon field. This can be eliminated, and the time derivative returned to its 
canonical form performing the Foldy-Wouthuysen transformation 

The illustration of the necessity for this field redefinition in the context of the heavy 
quark applications can be found in Ref. ||, Sect. 2.1. 

2.2 Applications to spectroscopy of heavy flavor hadrons 

To illustrate the consequences of the heavy quark Hamiltonian, let us consider the 
masses of hadrons containing a single heavy quark. This is not a dynamic question 
and requires only symmetry properties of Hq. It is described in much detail in a 
number of old reviews 0. 

The mass of a hadron Hq is given by the expectation value of the Hamiltonian: 



M h q = 7777 — { H Q\H tQ t\H Q ) 

and we have 



Ti-tot — flight + 7~Lq + rnq 

where we can expand 



with 



H Q = H + —Hx + \n 2 + ... (19) 
m Q m Q 



Hq = -ip^A ip Q ^ -Aj(0) [{flight}] 



—Hi = -—(7r 2 + aB) 



tuq 2m, 
1 „. 1 



T^ 2 — o 2 



(DE) + a ■ {E xn — nxE} , (20) 
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etc. Therefore, 



w x 1 (H Q \n 2 + aB\H Q ) - (j4-(i 2 g )h q 

MH Q = m Q + A + — ^- + - = m Q + A+ ^ + ... (21) 

The mass expansion was alternatively derived in Ref. [Q using the formalism based 
on the trace of the energy-momentum tensor. 

In Eq. (|2lD we introduced the notations fi 2 , /jg for the expectation values of two 
D = 5 heavy quark operators which will often appear in our discussion: 

= (H Q \n 2 \H Q ) QM = — — ^ (Hq\Qtt 2 Q(0)\Hq) qft 
zm h q 

fx 2 G = -(H Q \aB\H Q ) QM = ^l-{H Q \Q Z -a^Q(0)\H Q ) QFT . (22) 

The physical meaning of fi 2 is quite evident: the heavy quark inside Hq experiences 
a zitterbewegung due to its coupling to light cloud. Its average spatial momentum 
squared is fx 2 . The second expectation value measures the amount of the chromo- 
magnetic field produced by the light cloud at the position of the heavy quark. In 
principle, the actual heavy hadron states Hq depend on tuq ma the 1/mQ-suppressed 
terms of the Hamiltonian. Therefore, the above expectation values also have such 
terms. Often it is convenient to consider the asymptotic values at rriQ — > oo, and to 
use, correspondingly, the eigenstates of the mg^oo Hamiltonian TCq + flight- 

The parameter A appearing in Eq. (|21| ) was introduced as a constant in the Heavy 
Quark Effective Theory (HQET 0) in 0; it is associated with those terms in the 
effective Lagrangian £(//) (disregarded so far) which are entirely due to the light 
degrees of freedom. Needless to say that in the Wilsonian approach A is actually \x 
dependent, Wherever there is the menace of confusion the /x dependence of A 

will be indicated explicitly. 

There exists an interesting expression for this parameter through the anomaly 
in the trace of the energy- momentum tensor, derived in Ref. 

A = ^(HQ\^G 2 + J2rn q (l+ lm )q q \H Q ) mQ ->oo ■ (23) 

Here 7 m is the anomalous dimension of the light quark mass. This equation parallels 
the similar relation for the nucleon mass |7| well known in the chiral limit when all 
quark masses are neglected: 

M N = -±-(N\^G 2 \N) . 
2M N 4a s 

It should be noted, however, that the operator G 2 in Eq. (p3|), although local and 
gauge invariant, is not a local heavy-quark operator as they are understood in the 
framework of the heavy quark expansion: the latter must be of the form Q...Q 
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where ellipses denote a local operator at the same space-time point as the heavy 
quark fields Q. Quantum mechanically, the operators in Eq. (|23| ) describe the space 
integral over the whole volume occupied by the heavy hadron, rather than the value 
of the gluonic and quark fields at the position of the heavy quark. 

The renormalization properties of the operator G 2 are quite different in the 
sector of QCD with the heavy quark. An accurate consideration reveals [|]] that an 
additional term — //^p 2, must be added in the r.h.s. of Eq. (]23f) , where /i stands for 
the normalization point. The normalization point dependence of the heavy quark 
parameters will be discussed in more detail in the subsequent sections. 

The value of A = linim^oo ^M^-mQj has the scale of Aqcd and depends on 
the state of light degrees of freedom. These states generally carry spin j. In the 
limit vtlq — > oo the heavy quark spin decouples since the spin-dependent parts are 
present only starting 7i\. Thus, the heavy flavor hadrons can be classified not only 
by their total spin J but by the spin of light degrees of freedom j. It would be 
just the overall spin of a hadron in the hypothetical world with the spinless heavy 
quarks discussed in Sect. 1.1. Unless j = 0, there are two values of the total spin 
J = j ± |. The corresponding states form 'hyperfine' multiplets and are degenerate 
up to I/tuq corrections. They are, for example 

J = 1/2 | ^ * , j = 0,J=l/2 A c , A b 

For Ag-baryons all spin is carried by the heavy quark (up to I/tuq corrections). 
The observed spectroscopy of these states clearly supports this picture: 

M Ab - M B ~ 350 MeV M Ba - M B - ^ 90 MeV 

M Ac — Mb — 420 MeV M Ds - M D o ~ 104 MeV (24) 

These relations are easily improved including 1/toq terms, Eq. (0). The opera- 
tor Qtt 2 Q is spin-independent and its expectation value is the same for all members 
of a hyperfine multiplet. It does not split masses inside the multiplet. 

The chromomagnetic operator Q^aGQ = —2SqB(0) depends on the heavy 
quark spin Sq and thus lifts degeneracy leading to the hyperfine splitting among 
the members of the multiplet. Restricted to a particular hyperfine multiplet, the 
chromomagnetic field 5(0) is proportional to the spin of light degrees of freedom: 
B = c ■ j. Therefore, 

(ffB) = 2c(S Q j) = c(j(J+l)-jU + l)-l) ■ (25) 
It is easy to see that 

£ (H Q \QaBQ\H Q ) = Tr QaBQ = (26) 

always holds if the summation is performed over a hyperfine multiplet. Therefore, 
for example, 

+ 3/4(5*) = 0. (27) 
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In the A;, family the expectation value of b^aGb vanishes. 

So far most of the practical applications refer to B mesons; as a result, usually 
/•4 P ro P er denotes the expectation value of the kinetic operator just in B or B*. 
Likewise, [Mq = fJ-a(B) = —3[1q(B*) . Experimentally, 

M B ,-M B ~ (I + l) A = !A^46MeV. (28) 
\3 / 2m b 3 2m b 

Neglecting the difference between M B + M B * and 2m b (which is formally an effect 
of higher order in l/m b ), one can write 

ill ~ ^(M B ,-M B ) ~ 0.36 GeV 2 . (29) 

In charmed mesons Md*—Md ~ 140 MeV which agrees with the fact that this hyper- 
fine splitting is proportional to l/m c - The hadron mass averaged over a hyperfine 
multiplet, e.g. M B = 3Mb *+ Mb - IS a ff ec ted at order 1/uiq by only the kinetic energy 
term /^/2m&. 

The mass expansion can be extended to higher orders in I/tuq. For example, to 
order I/tuq the hadron mass takes the form || 



mq=oo 



(30) 

The nonlocal (positive) correlator p 3 is the second-order iteration of the 1 / mq part 
of the Hamiltonian: 

P 3 = *J d 4 x -^-(H Q \T{Q{airfQ{x), Q(^) 2 Q(0)}\H Q y mQ=OD (31) 

(the prime indicates that the diagonal transitions must be removed from the cor- 
relation function). In B mesons the expectation values of the two terms in l/rng 
Hamiltonian, the Darwin and the convection current interactions, are denoted as 

As = ^(B\b3-ExHb\B) = -JL-(B'\b<j-Ex9b\B') . (32) 

The Darwin operator by virtue of QCD equations of motions equals to the local 
four-fermion operator of the form —^fQ^QY^q^lQ^Q- 

A similar decomposition can be made for the nonlocal correlators whose expecta- 
tion values were generically denoted by p 3 . In particular, in pseudoscalar and vector 
ground states 

fL = if d 4 x -^(B\T{W%{xl H 2 b(0)}\BY 
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Pla = if d'x ^L-(B\T{bn 2 b(x), baBb(0)}\B)> 

^ijSkt + ^ A (S ik 6 jl -6 u 6 jh )=i J d 4 x j^{B\T{baiB k b{x), ba^b^B)' . 

(33) 

Then one has for p 3 in B and B*, respectively, 

(P 3 ) B = PL + Pla + Ps + Pa , (P 3 ) B , = PL - \pIg + P%- \p\ ■ (34) 

These parameters also give the 1/rriQ corrections to the expectation values of the 
kinetic and chromomagnetic operators Q. 

In many applications one needs to know the difference between m b and m c . This 
difference is well constrained in the heavy quark expansion. For example, 

M B + 3M B * M D + 3M D , u 2 /l 1\ P 3 D -p 3 (l l\ /l' 
m b -m c = + — + — - — — ^ + ° — 3 

(35) 

Here is the asymptotic expectation value of the kinetic operator and p 3 = p 3 7r + 
p| is the sum of two positive nonlocal correlators. As will be discussed in the 
subsequent sections, all quantities in Eq. (|35|) but the meson masses depend on the 
normalization point which can be arbitrary, except that it must be much lower than 
m Cj fe. Evaluating the above expansion we arrive at 

u 2 _n 5 GeV 2 

m b -m c ~ 3.50 GeV + 40 MeV ^ — — + AM 2 , I AM 2 | £ 0.02 GeV, (36) 

0.1 GeV 2 

where plausible assumptions about the D = 6 expectation values have been made. 
The values of the hadronic parameters will be discussed later. Let us mention that 
the m b —m c estimate at p 2 ~ 0.5 GeV 2 appears to be in a good agreement with the 
separate determinations of m b and m c from the sum rules in charmonia and T. 



2.3 Heavy quark symmetry for formfactors 

The amplitudes of the semileptonic weak b —>■ c decays are described by the corre- 
sponding transition formfactors. Typical semileptonic b — >• c decays as they look like 
in Feynman diagrams are shown in Figs. 1. The hadronic part of the weak decay 
Hamiltonian mediating such decays is 

H weak = J d 3 xe-^c llM (l- l5 )b(x) . (37) 

It says that b with a momentum p is instantaneously replaced by the c quark with the 
momentum p — q. The resulting state hadronizes into eigenstates of the Hamiltonian 
corresponding to mg=m c , that is, must be projected onto such states. 

The space-time picture of the decay is simple for heavy quarks. At t < the 
initial b hadron is at rest and constitutes a coherent state of light degrees of freedom 
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Figure 1: Quark diagrams for the exclusive B — > £)(*) (a) and generic (b) semileptonic 
decays. 

in the static field of the heavy quark. At t — the b quark emits the lepton pair 
with the momentum q and transforms into a c quark. The c quark gets the recoil 
momentum — q and starts moving with the velocity v = —q/m c . Such a state is 
not anymore an eigenstate of the Hamiltonian, and afterwards undergoes nontrivial 
evolution. The light cloud can get a coherent boost along the direction of — q and 
form again the same ground-state, or excited hadron. Alternatively, it can crack 
apart and produce a few-body final hadronic state. 

The heavy quark symmetry per se cannot help calculating the amplitudes to 
create such final states. However, it tells one that the hadronization process does 
not depend on the heavy quark spin, or on the concrete value of the mass uiq but 
rather on the velocity of the final state heavy hadron. This velocity cannot change 
in the process of hadronization if only soft gluons are exchanged between the heavy 
quark and the light cloud. This independence holds, of course, only when uiq is 
very large (and the final state quark does not move too fast); there are various 
1/rriQ corrections at finite masses. 

Let us consider, for example, the ground state transition B — > D. Its amplitude 
depends on the velocity v of D, f{v 2 ). The very same function would describe also 
decays B — > D*, or the elastic amplitude of scattering of a photon on the b quark 
B — > B, B — > B* . Moreover, in the proper normalization /(0) = 1 holds. This fact 
follows from the conservation of the 6-quark vector current, for the amplitude at 
zero momentum transfer measures the total '6-quark charge' of the hadron n^ — ni. 
Its origin is simply understood: if {7=0, the final state is not really disturbed by 
movement of the static source. A nontrivial rearrangement of the light cloud for 
B — > D is associated only with the mass-dependent terms in Hq vanishing when 
rriQ -> oo. 

Consider the vector bq^b current in B meson. It is described by the single 
formfactor f+(q 2 ): 

(B{p')\hm\B{p)) = U(q 2 )(P + P'),, %=V»-V'» (38) 

(the second structure is forbidden by T invariance or current conservation). 

The value /+(0) measures the total 'beauty charge' of the hadron and is not renor- 
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malized by the strong interaction, /+(0) = 1.0 Passing to the velocities, we use 



instead of q 2 the scalar product Vnv' : 



yuv ) = — 

and 



™') = ^=1-^>1, (39) 



M| 2M B 



f + (q 2 ) = £(«/) , e(l) = 1 • (40) 
£ (zw') is called the Isgur-Wise function. The heavy quark symmetry then states that 

w)|^<°)i*m> - (^§gfe<*+^ - imm**-^ « W) ■ (41) 

.D* differs from .D only by the alignment of the c quark spin. Taking this into account 
yields 

(D*(v',e)\H^b(0)\B(v)) = ze^eX^K) ^M B M* D 

(D*(v',e)\c W b(0)\B(v)) = fe(vv'+l) - v'^v)} t{vv') sjM B M* D . (42) 

It is important that these relations are valid in the limit m^ c — > oo and if no short- 
distance radiative corrections were present. The corrections to the symmetry limit 
are minimal at q — (v — v', the so-called zero recoil point); numerically the value 
of the axial formfactor was estimated to be Fb^d*(0) ~ 0.9 @. The corrections at 
arbitrary v ~ 1 are generally significant. 

2.4 Feynman rules at tuq —> oo 

The Feynman rules for heavy quarks are usual propagators and vertices for nonrela- 
tivistic particles where 1/m — *■ 0. If the heavy quark momentum is p^ = (rriQ+u), p) 
then 

G(p) = lim 5 a p ^-Sij = — — 

2m 

= g.tf-6,a . (43) 

Here a, a, (3 are color indices and spinor indices i,j take values 1 or 2. (It is 
often advantageous to keep the nonrelativistic term p 2 /2rriQ in the propagator as 
an infrared regulator.) These rules follow immediately from the static Lagrangian 



a 



V 

£q = PQiDaVQ , D = d - iAo— . 

The same nonrelativistic system can be considered in an arbitrary moving frame, 
where one can write 

C v = ifti{vD)ip v . (44) 

4 For simplicity, we use the convention that B consists of b and q. That is, B~ is a B meson 
while B + is B. 
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Instead of the nonrelativistic spinor ipq one then considers the "bispinor" <p v (x) = 
■^■Q(x); the propagator is written as 

1+j m^^-^ 1+^ _^ 1+^ 1 
2 m 2 Q -(p-k) 2 -ie 2 2 wfc-ze 1 ' 

y 

where fc = p—rriQV. Of course, such a generalization can be useful only if the initial 
and final state hadrons have different velocities, v ^ v'. In that case the external 
('weak') source carries a large momentum q ~ rriQSv. 

2.4.1 Subtleties of actual QCD 

It must be noted, however, that the field-theoretic description of processes where the 
heavy quarks change velocity becomes subtle if the quantum radiative corrections 
are really incorporated. Any change in the velocity of a static source leads to actual 
radiation of real hadrons with momenta k all the way up to tuq: 



dw Oi s {uj) , 2 



duj uj 



(v'-vf ; (46) 



here uj denotes the radiated energy. This means that with infinitely heavy quarks 
one would actually observe radiation of light hadrons off the heavy flavor hadrons 
with arbitrary large energies. Without an ultraviolet cutoff any exclusive transition 
probability would vanish being suppressed by the more or less universal factor (the 
square of the non-Abelian analogue of the nonrelativistic Sudakov formfactor) 

S ~ e -1^) 2ln 7 . (47) 

Here /i denotes the ultraviolet cutoff and e determines the energy resolution for the 
final hadronic state (for the validity of the perturbative expansion e must be much 
large than Aqcd; taking it around the typical hadronic scale gives the estimate of 
the overall perturbative suppression of the exclusive transition probability). This 
factor has a meaning of the probability of the heavy color source not to emit gluons 
in the interval of energies between e and \x in the act of acceleration. 

Beyond the small velocity approximation the factor (Av) 2 in Eq. ( |4"7[ ) must be re- 
placed by | (t^=t In ^zj^j — 2^ which can be read off the expression for the intensity 
of electromagnetic radiation in classical electrodynamics: 



1 dl(u) 
uj du 




(48) 



It is remarkable that this law is not modified by quantum corrections f| as long as 
uj remains much smaller than the masses of all charged particles and the effects of 
their virtual production can be neglected. 
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In the non-Abelian theories like QCD the situation is different, even in the dipole 
approximation \v\ <C 1. The non-Abelian QCD dipole radiation was considered in 
detail in Ref. [fTOfl . In contrast with QED where the coupling for the soft photon 
radiation is exactly a em (0) and is not renormalized by quantum corrections, in QCD 
the dipole radiation is governed by the effective coupling a^(to) which runs with 
the energy scale: 

with a s the MS coupling. It can be shown that the second (conformal) term in the 
O(of) part of the coupling must equal twice the corresponding part of the so-called 
cusp anomalous dimension of the Wilson lines (at small angle) investigated in detail 



in the mid 80's [fLi]] . The dipole radiation coupling governs the normalization point 
evolution of the heavy quark masses and a number of heavy quark operators in the 
effective low-energy theory 

As any effective coupling, q^w) is not a purely perturbative object at arbitrary 
precision. It has a power-suppressed nonperturbative component at arbitrarily large 
energy u which, for example, depends on the particular type of the heavy flavor 
hadron which is accelerated. (Pure perturbation theory does not depend on the 
light cloud surrounding heavy quark for actual hadrons.) Using the OPE approach 



it was shown [10| that such nonperturbative effects in the dipole radiation fade out 



at least as 1/cu 3 , and for the ground state hadrons (B mesons) it was estimated that 

KV))^-^) 3 - ™ 

In general, the perturbative factors suppressing the velocity-changing transition 
amplitudes of infinitely heavy quarks are universal and related to the cusp anomalous 
dimensions of Wilson lines describing propagation of massive color objects. They 
appear in the renormalization of bent Wilson lines; the cusps correspond to instant 
changes in velocity and thus lead to new ultraviolet divergences. The dedicated dis- 
cussion can be found in the original publications RTI| . The nonperturbative aspects 



of factorization in the radiative effects have not been carefully studied yet. 



3 Basic Parameters of the Heavy Quark Expan- 
sion 

3.1 The heavy quark mass 
3.1.1 What is m Q l 

In quantum field theory the object we begin our work with is the Lagrangian for- 
mulated at some high scale Mq. The mass mo is a parameter in this Lagrangian; 
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it enters on the same footing; the bare coupling constant with the only 

difference being that it carries dimension. As with any other coupling, m enters all 
observable quantities in a certain combination with the ultraviolet cutoff M , which 
is universal for a renormalizable theory. 

The mass parameter mo by itself is not observable, like . For calculating 
observable quantities at the scale /i <C M h is usually convenient to relate m to 
some mass parameter relevant to the scale fi. Integrating out momentum scales 
above /i converts into a s (fx) - and likewise m into m{jj). Such m(/z) is not 
something absolute since depends on \i. It is either used on the same footing as a s (/i) 
or, in the final expressions, is eliminated in favor of some suitable observable mass. 
For example, in quantum electrodynamics (QED) at low energies (i.e. E <C m e ) 
there is an obvious "best" candidate: the actual mass of an isolated electron, m e . 
In the perturbative calculations it is determined as the position of the pole in the 
electron Green function (more exactly, the beginning of the cut). The advantages 
are evident: m e is gauge- invariant and experimentally measurable. 

The analogous parameter for heavy quarks in QCD is referred to as the pole 
quark mass, the position of the pole of the quark Green function. Like m e it is gauge 
invariant. Unlike QED, however, the quarks do not exist as isolated objects (there 
are no states with the quark quantum numbers in the observable spectrum, and the 
quark Green function beyond a given order has neither a pole nor a cut). Hence, 
m pole cannot be directly measured; m pole exists only as a theoretical construction. 

In principle, there is nothing wrong with using m polc in perturbation theory where 
it naturally appears in the Feynman graphs for the quark Green functions, scattering 
amplitudes and so on. It may or may not be convenient, depending on concrete goals. 

The pole mass in QCD is perturbatively infrared stable, order by order, like in 
QED (the formal proof was given recently in [[HJ). It is well-defined to every given 
order in perturbation theory. One cannot define it to all orders, however; the sum 
of the series does not converge to a definite number. In a sense, the pole mass is 
not infrared-stable nonperturbatively. Intuitively this is clear: since the quarks are 
confined in the full theory, the best one can do is to define the would-be pole position 
with an intrinsic uncertainty of order Aqcd |~3| . 



Based on experience in QED or ordinary QM, non-existence of mg lc may seem 
counter-intuitive. Employing perturbation theory, we start with the free quark prop- 
agator 

G{p) ' 



(mg(/i) is the parameter entering the Lagrangian) which has a pole at p 2 = mg(/j) 
and, therefore, describes a particle with mass mq{p:). Accounting for the gluon 
exchanges to the first order in a s adds the diagram Fig. 2 and the pole moves to 

p 2 ~ (niQ(fi) + ^-f^j describing now a particle with the mass itlq ~ mg(/i) + ^// , 
etc. In any order of perturbation theory we see a quark pole and the corresponding 
particle with certain mass, differing from the mass in the Lagrangian. This mass is 
just the "pole" mass. Clearly, it depends on the order of perturbation theory one 
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considers, and on a concrete version of the employed expansion: 
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(51) 



n=0 



It is tempting to define the 'actual' pole mass as a sum of the series fl5"l"|). It 
appears, however, that the sum does not converge to a reasonable number and 
cannot be defined with the necessary accuracy in a motivated way, but suffers from 
an irreducible uncertainty of order Aq CD . 

Before explaining the origin of this perturbative uncertainty, let us remark that 
such a definition of the quark mass is nothing but equating it with the lowest eigen- 
value of the (perturbative) Hamiltonian in the sector with the number of heavy 
quarks 1. The eigenvalues of the Hamiltonian, however, are not short-distance quan- 
tities, and their calculation in the perturbative expansion is not adequate. This is in 
the nice correspondence with the observed fact that no isolated heavy quark exists 
in the spectrum of hadrons which contains instead B, B**, B s . . .with masses 
differing by amount C(Aq CD ). 

The physical reason behind the perturbative instability of the long-distance 
regime is the growing of the interaction strength a s . One can illustrate this in- 
stability in the following way. Consider the energy stored in the chromoelectric field 
in a sphere of radius R 3> I/tuq around a static color source, 



This energy is what one adds to the bare mass of a heavy particle to determine 
what will be the mass including QCD interactions. The definition of the pole mass 
amounts to setting R — > oo; i.e., in evaluating the pole mass one undertakes to 
integrate the energy density associated with the color source over all space assuming 
that it has the Coulomb form. In real life the color interaction becomes strong at 
Rq ~ 1/Aqcd; at such distances the chromoelectric field has nothing to do with the 
Coulomb tail. Thus, one cannot include the region beyond Rq in a meaningful way. 
Its contribution which is of order Aq C d, thus, has to be considered as an irreducible 
uncertainty which is power-suppressed relative to tuq, 



Exactly this behavior is traced formally in perturbation theory. In the nonrel- 
ativistic regime where the internal momentum \k\ <C tuq the expression for the 
diagram Fig. 2 is simple, 




(52) 




(53) 




(54) 
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Figure 2: Perturbative diagrams leading to renormalization of the heavy quark mass. 
The contribution of the gluon momenta below tbq expresses the classical Coulomb self- 
energy of the colored particle. The number of bubble insertions into the gluon propagator 
can be arbitrary generating corrections in all orders in a s . The factorial growth of the 
coefficients produces the IR renormalon uncertainty in m^ c of order Aqcd- 



The latter expression is |V(0) with V(R) the usual Coulomb potential between two 
quarks. The running of the coupling is generated by dressing the gluon propagator 
by virtual pairs and leads to 

4 rd z ka s {P) 
6m ^3J^—- (55) 

Since 

« s (fc 2 ) = a s ( / i 2 )|l + ^^Mn^| , b=jN e -ln f , (56) 

we can expand a s {k 2 ) in a power series in a s (fi 2 ) and easily find the (n+l)-th order 
contribution to Smq , 

H n+1) „t?M n \ ( b <dY , (57) 

m Q 3 7r \ 2n J 

The coefficients grow factorially and contribute with the same sign. Therefore, 
one cannot define the sum of these contributions even using the trick with the Borel 
transformation. The best one can do is to truncate the series judiciously. An optimal 
truncation leaves us with an irreducible uncertainty ~ (9(Aq CD ) |i~5| . The above 
perturbative corrections are example of the so-called infrared renormalons [ R], [T7J . 



This uncertainty can be quantified. A formal Borel resummation of such non- 
summable series leads to the result which literally has an imaginary part which can 
be taken as a measure of the uncertainty. The imaginary part for the series Eq. 
is 

3tt 

It became conventional to assign the formal imaginary part divided by it to the 
irreducible uncertainty. Even with this minimal choice 



Im< = ^e 5 / 6 AgS D . (5* 



<r pole 1 
^ 7T 



Immf = ^e 5 / 6 A^ D - 0.7A^ D . (59) 



27 
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Thus, the perturbative expansion per se anticipates the onset of the nonpertur- 
bative regime (the impossibility of locating the would-be quark pole to accuracy 
better than Aqcd)- Certainly, the concrete numerical value of the uncertainty in 
m poie obtained through renormalons is not trustworthy. The renormalons do not 
represent the dominant component of the infrared dynamics. However, they are a 
clear indicator of the presence of the power-suppressed nonperturbative effects, or 
infrared instability of m pole ; the very fact that there is a correction (9(Aqcd/^q) is 
beyond any doubt. 

It is worth noting that the pole mass was the first example where a quantity which 
is perturbatively infrared- stable was shown not to be stable nonperturbatively at the 
level Aq CD . The observation of Refs. [[14], [15| gave impetus to dedicated analyses of 
other perturbatively infrared- stable observables in numerous hard processes without 
OPE, in particular, in jet physics. Such nonperturbative infrared contributions, 
linear in Aqcd/Q were indeed found shortly after in thrust and many other jet 
characteristics (for a review and a representative list of references see e.g. HTTP). 

Since it is impossible to relate mQ(fi) and TriQ' le to the necessary accuracy, it 
is clear that either mg le or mg(/i) must be irrelevant for the 1/m.Q expansion, 
for example, for calculating the decay widths of heavy flavors which depend on 
high powers of the quark masses. Which mass must be used then? This question 



was formulated and answered in the early 1994 in Ref. [14|. In agreement with 
the qualitative discussion given above, the answer is: rn^ e is irrelevant and must 
be replaced everywhere by rriQ^fi). For the inclusive decay widths, this can be 
understood as the perturbative counterpart of the QCD theorem [|Tj| stating that 
the mass gap ~ C(Aq CD ) differentiating the hadron mass from the quark mass, does 
not affect the width - the pole mass dressed by soft interactions is the perturbative 
analogue of the hadron mass. The 1 /irtq infrared renormalon uncertainty in the 
pole mass disappears from the perturbative corrections when the width is expressed 
in terms of the short-distance heavy quark mass; it is present when one attempts to 
compute the width through the pole mass. 

Another transparent way to illustrate irrelevance of the pole mass is to vary the 
number of space dimensions, viz. to descend to the (2+1) theory. In D = 3 the 
(pole) mass logarithmically diverges in the infrared, cf. Eq. (0). The OPE, at the 
same time guarantees that the total widths are infrared safe in any dimension, so 
this infrared part cannot affect the width. The infrared divergence of the mass in 
D = 3 is simply reflection of the Coulomb potential In R growing at large R. It is 
evident that total decay probabilities of a heavy quark must know nothing about 
the behavior of the interaction at infinite distance. 

In the OPE, the infrared part of the pole mass is not related to any local operator 
QOiQ, and does not enter any observable calculable in the short-distance expansion. 
Since this infrared piece does not enter observables, it cannot, in turn, be determined 
experimentally. The numerical instability of various attempts to pinpoint the value 
of ?T7<Q 0le is a result of using the perturbative expansion for the effects originating from 
the nonperturbative domain. The above facts were later illustrated in Ref. |lll in a 



concrete model for the higher-order corrections to the semileptonic widths obtained 
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via the so-called bubble resummation of the one-loop perturbative diagrams. 

It is important to note that the irrelevance of the pole mass goes beyond the 
problems with the infrared renormalon contributions illustrated above. Even if there 
existed some way to define reasonably the sum of the pure perturbation series for 
mg° le , or asymptotic states with single-quark quantum numbers with finite energy 
existed in a strong-interaction theory like QCD, this mass still would have been 
inadequate for constructing the effective field theory, to the extent the difference 
with mg(/i) cannot be neglected numerically. 

Summarizing, the pole mass per se does not appear in OPE for infrared-stable 
quantities. Such expansions operate with the short- distance (running) mass. Any 
attempt to express the OPE-based results in terms of the pole mass creates a prob- 
lem making the Wilson coefficients ill-defined theoretically and poorly convergent 
numerically. 

A properly constructed perturbative treatment suitable for the 1 / tuq expansion 
incorporates only gluons with \k\ J> /x which are not 'resolved' and are included into 
the heavy quark field wavefunction corresponding to the heavy quark field Qq^(x) 
normalized at the scale /x. The normalization point /x can be changed: descending 
from /x down to /xi < /x one has to integrate out newly-unresolved gluons with 
\L\ < \ k\ < /x. For example, the Coulomb field associated with such gluons increases 
the mass of the quark by the amount 



*»<.= / ,. («) 



d 3 k 4 a s (k 2 ^ 
4tt 2 3 p 

The pole mass clearly appears when /xi — > 0. From the OPE point of view, it 
is an attempt to construct an effective theory with the normalization scale /x = 
formulated, nevertheless, still in terms of quarks and gluons. Speaking theoretically, 
one can imagine a limit of small /x which would correspond to integrating out all 
modes down to /x = in evaluation of the effective Lagrangian. It would be nothing 
but constructing the S'-matrix of the theory from which one can directly read off all 
conceivable amplitudes. Clearly, it could have been formulated in terms of physical 
mesons and baryons but not of quarks and gluons. 

In the Wilson OPE one uses mg(/x) with Aqcd <C /x <C m&. We illustrate later 
that just such a mass can be accurately measured in experiment. It is /x-dependent: 

dra Q (/x) 16 a 8 (n) 4 a,(/x) /x f 2 fi 2 \ 

— i = - IT o + ° \ a si a s— ; (61) 

d/x 9 7r 6 ix rag \ tUq J 

the higher order perturbative corrections were computed recently [[ll]]. There are 
different schemes for defining mg(/i) (similar to renormalization schemes for a s ), and 
the coefficients above are generally different there. As long as a concrete scheme 
is adopted, there is no ambiguity in the numerical value of ttt.q(/x). Instead of the 
HQET parameter A in QCD one has A(/x) = lim^^oo Mn Q — mg(//). The value of 
A(/x) is of the hadronic mass scale if /x does not scale with uiq. 

There exists a popular choice of a short- distance mass, the so-called MS mass 
m(/x). The MS mass is not a parameter in the effective Lagrangian; rather it is a 
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certain ad hoc combination of the parameters which is particularly convenient in 
the perturbative calculations using dimensional regularization. Its relation to the 
perturbative pole mass since recently is known already to three loops fl20f : 



m? e = m Q (m Q ) (l + ^^^ + (1.566- 3.73) f^V + ....) (62) 



3 71 \ 7V 

At \i > rriQ the MS mass coincides, roughly speaking, with the running Lagrangian 
mass seen at the scale ~ /i. However, it becomes rather meaningless at fi <C vtlq\ 



m Q {fi) ~ m Q (m Q )\l-\ In \ . (63) 




It logarithmically diverges when fi/mq — > 0. For this reason rh(fi) is not appro- 
priate in the heavy quark theory where the possibility of evolving down to a low 
normalization point, fi <C rriQ, is crucial. Otherwise, for example, Mu Q —fnQ oc rriQ 
and does not stay constant in the heavy quark limit. 

The reason for this IR divergence is that the MS scheme technically attempts to 
determine the (perturbative) running of all quantities by their divergences calculated 
when the space-time dimension approaches D — 4. The divergence can emerge 
only at k — > or k —>■ oo. For the mass, the IR divergences are absent, and the 
dimensional regularization is sensitive to the UV divergence at k 3> rriQ. Studying 
only 1/(D — 4) singularities, it is unable to capture the change of the regime at 
fi ^ vhq and assumes the same running in this domain. The actual running below 
rtiQ is slower. 

The properly defined short- distance masses always exhibit an explicit linear /i- 
dependence similar to Eq. ( |B"TD at fi <C uiq. The perturbative pole mass, order by 
order, would correspond to the limit fi — > 0. However such a limit does not exist. 

Since rn V Q le does not exist as a well-defined mass parameter, a different, short- 
distance mass must be used. The normalization point fi can be arbitrary as long 
as fi ^> Aqcd- It does not mean, however, that all masses are equally practical, 
since the perturbative series are necessarily truncated after a few first terms. Using 
an inappropriate scale makes numerical approximations bad. In particular, relying 
on rfiQ(m,Q) in treating the low-scale observables can be awkward. The pedagogical 
example illustrating this point can be found in Refs. El], Ej . 



Needless to say, it is the high-scale masses that appear directly in the processes 
at high energies. In particular, the inclusive width Z — > bb is sensitive to m^iMz)] 
using MS mass normalized at fi ~ Mz is appropriate here. On the contrary, the 
inclusive semileptonic decays b — > civ are rather low-energy in this respect WA \ this 
is true even for b —>■ u. 

The construction of the running low-scale heavy quark mass suitable for the OPE 
to any order in perturbation theory can be done in a straightforward way. However, 
the Wilsonian approach implies introducing a cutoff on the momenta of the gluon 
fields. Since gluon carries color, its momentum is not a gauge-invariant quantity, 
and such a mass typically looks not gauge-invariant. More accurately, obtaining the 
same m^/i) requires somewhat different cutoff rules in different gauges. 
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Even though this is not a real problem for the theory, it is often viewed as a 
disadvantage. To get rid of this spurious problem, a manifestly gauge-invariant 
definition of the running mass was suggested in |^2, ^| which is formulated only 



in terms of observables. The idea is to explicitly subtract from the pole mass the 
infrared pieces it contains. These are uniqually determined by the so-called small 
velocity (SV) sum rules in the heavy quark limit discussed in Sect. 4. They are 
moments of the structure functions of the infinitely heavy quark in the process 
where its velocity changes by a small amount. In the presence of hard gluons these 
moments diverge in the ultraviolet, and must be cut off at some energy /i. This 
cutoff enters then as the normalization point for the heavy quark mass. 

Such defined mass is convenient for the OPE in the 1 /mq expansion since more 
or less directly enters many relevant processes, e.g. heavy flavor transitions or the 
threshold heavy flavor production. Its relation to the MS mass is known with enough 
accuracy JTOf- To avoid ambiguities we always use this definition unless other con- 
vention is indicated explicitly. 

It is important to note the following fact. In the relativistic theory for a particle 
with mass m one always has p 2 = m 2 , that is, E — ^m' 2 +p 2 . In the nonrelativistic 
expansion, therefore, 

E = m+£--£ + ... ; (64) 

all coefficient are fixed in terms of powers of m by Lorentz invariance. In applications 
to heavy quarks it is advantageous to use such a Wilsonian cutoff which preserves 
usual QM properties for the price of apparently violating Lorentz invariance. This is 
justified since the problem from the very beginning has physically preferred frame - 
the one in which the heavy quark is at rest. (More detailed discussion can be found 
in |p3|.) Then the mass parameters in Eq. ([64]) generally become different: 



H Q = mo(p%(p Q - (p^A (p Q + — (iiDf + c G aB^j - ... (65) 

even for a "quasifree" quark in the effective theory. This phenomenon is well known 
in the solid state physics. Of course, the difference between m, appears only due to 
perturbative corrections: 

mi(ii) - m k (ii) = 0{a s n) . 

Moreover, the differences can be calculated perturbatively and are completely free 
from any infrared effects existing below the scale /i. For example, in our scheme 
m (/i) — 777.2 (/•*) + | ■ We use the mass 777,2 (/i), the mass that enters the kinetic 
energy operator p 2 /2m; it is the most relevant mass for heavy quark transitions. 

This calculable difference between different "masses" for the same quark must be 
properly accounted for in the OPE analysis. Probably, the most obvious place where 
it plays a role is the QQ threshold physics. While in the short- distance expansion 
the free quark threshold starts at 2mo(ju), the propagation of heavy quarks or the 
bound state dynamics are actually determined by m 2 (/u). The shift in the position 
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of the threshold which serves as a reference point for energy in the nonrelativistic 
system must be properly taken into account. 

Concluding this section, let us make a side remark concerning the t quark mass 



2jJ . The peculiarity of the t quark is that it has a significant width T t ~ lGeV 
due to its weak decay. The perturbative position of the pole in the propagator 
is, thus, shifted into the complex plane by — |T t . The finite decay width of the 



t quark introduces a physical infrared cutoff for the infrared QCD effects [24]. In 
particular, the observable decay characteristics do not have ambiguity associated 
with the uncertainty in the pole mass discussed above. The uncertainty cancels in 
any physical quantity that can be measured. That is not the case, however, in the 
position of the pole of the t-quark propagator in the complex plane (more exactly, 
its real part). The quark Green function is not observable there, and one would 
encounter the very same infrared problem and the same infrared renormalon. The 
latter does not depend on the absolute value of the quark mass (and whether it is 
real or have an imaginary part). Thus, in the case of top, one would observe an 
intrinsic - but artificial - infrared renormalon uncertainty of several hundred MeV 
in attempts to relate the peak in the physical decay distributions to the position of 
the propagator singularity in the complex plane. 

3.1.2 The numerical values of m c and m;, 

The mass of the c quark at the scale ~ m c ~ 1 GeV can be obtained from the 
charmonium sum rules P5|1 , rh c (m c ) ~ 1.25 GeV. The result to some extent is 
affected by the value of the gluon condensate. To be safe, we conservatively ascribe 
a rather large uncertainty, 

m c (m c ) = 1.25 ±0.1 GeV. 

There are reasons to believe that the precision charmonium sum rules actually de- 
termine the charmed quark mass to a better accuracy. 

An accurate measurement of m& is possible in the bb production in e + e~ anni- 
hilation. Since we want to know mj, with comparable or better absolute precision, 
both the data and calculations, at first glance, must have increased accuracy. The 
data are available, however, only below and near the threshold. Certain integrals 
(moments) of the cross section over this domain are particularly sensitive to the 
low-scale mass m^fi) with /x in the interval 1 to 2 GeV. 

Using dispersion relations 



iW) = n s( o) + Jl / ™ (66) 

one evaluates the polarization operator Hb(q 2 ) (and its derivatives) induced by the 
vector currents by^b, in the complex q 2 plane at an adjustable distance A from the 
threshold. Such quantities are proportional to weighted integrals over the experi- 
mental cross section; the integrals are saturated in the interval ~ A near threshold, 
and are very sensitive to the mass m^A). 
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A dedicated analysis of this type was first carried out by Voloshin who con- 
sidered a set of relatively high derivatives of Uf, at q 2 = 0. On the phenomenological 
side they are expressed through moments of Rb(s), 




(67) 



while the theoretical expressions for the very same moments are given in terms 
of the b quark mass and a s . The relevant momentum scale here is « ~ m&/ ' \Jn. 
Considering small-n moments I n one would determine m& at the scale of the order 
m&. The small-n moments are contaminated by the contribution of Rb above the 
open beauty threshold where experimental data are poor. Increasing n shrinks the 
interval of saturation and, thus, lowers the effective scale. On the other hand, 
we cannot go to too high values of n where infrared effects (given, first, by the 
gluon condensate) explode. There is still enough room to keep the gluon condensate 
small and, simultaneously, suppress the domain above the open beauty threshold. 
In the fiducial window, u must be large enough to ensure control over the QCD 
corrections. The latter requires a nontrivial summation of enhanced Coulomb terms 
unavoidable in the nonrelativistic situation. As known from textbooks, the part 
of the perturbative corrections to the polarization operator, associated with the 
potential interaction, is governed by the parameter a/\v | rather than by a per se. 

Let us briefly illustrate why the moments I n determined in this way pinpoint the 
running mass r?v Naively, the cut in Ui,(q 2 ) starts in the perturbative calculations 
at 2m^ ole , which seems to determine the strongest dependence on the b quark mass. 
However, the bb production is also affected by the potential interaction of b and b; the 
latter even generates a number of bound states below the bb threshold. The gluon 
exchanges increas m^ ole and thus tend to suppress /„ (we keep the short distance 
mass fixed). However, the Coulomb effects enhance the moments both owing to the 
emerging bound states and due to attraction above the threshold. 

The largest infrared contribution to the pole mass is linear in the gluon momen- 
tum, Eq. (0). Its cancellation can be understood on the example of QED with 
the Abelian gauge interaction. The expression for the mass shift is simply self- 
interaction |Vir(0) where Vjr is the heavy quark potential mediated by the gauge 
interactions with momenta below certain u <C nib HU, [14]]. Yet the mass of the bb 
system includes also the same Coulomb interaction between quark and antiquark. 
Since for the colorless bb state the sum of color "charges" is zero, these effects can- 
cel each other for the quanta with wavelength less than the interquark spacing r. 
Therefore, for the Fourier transform of the potential V(q ) in terms of which 



only the components with \q \ > 1/r contribute. The softer exchanges are suppressed 
by powers of the multipole factor q 2 r 2 . 




(68) 
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This, of course, automatically emerges in all calculations. Let us single out the 
effect of gluon exchanges with \q\ < jj,: 

= - L„ S v ^ e " p = - y ° + r v v '-- (e9) 

v =f £%v(i), v 2 =f ^-Xm), ... 

(the minus sign reflects the fact that the second particle is an antiquark). If quarks 
reside at distances much smaller than the soft potential is just a constant. Its 
sole role is only to shift the energy of all bb states by a constant amount — Vq. It 
does not affect wavefunctions and, therefore, does not modify the coupling of the 
virtual photon in e + e~ annihilation to these states. 

However, a constant potential A cannot change the energy of the neutral system, 
whether the field is classical or quantum. It means that just the opposite shift is 
present in the sum of masses of b and b renormalized by the same gauge interaction, 
which is self-manifest in Eq. d54|) . The above arguments show that the cancellation 
of the infrared effects in the quark mass is a nontrivial consequence of the gauge 
nature of QCD interactions, the fact explicit in the OPE. 

The above reasoning, while illustrating the cancellation of the infrared effects 
present in the pole mass, is too simplified in many aspects. A more dedicated 
discussion can be found in Ref. M. 

Over the last few years the perturbative decsription of the moments I n has been 
significantly improved and now includes the next-to-next-to-leading (NNLO) terms. 
The corresponding analysis of the data was performed by a few groups [p^] , differing 
in technical details; the values of m^l GeV) were found around 4.57 GeV with the 
estimated uncertainty 40 to 60 MeV. In what follows we accept the value of the 
running b quark mass 

m b (l GeV) = 4.57 ± 0.05 GeV (70) 



3.2 fil and /Iq 

The heavy quark masses m c , mi,, being the key parameters in the HQE, are to a large 
extent 'external' to the properties of the effective low energy theory itself. There 
are two nonrelativistic heavy quark operators in the Hamiltonian; their expectation 
values Eqs. (|22|) in the heavy meson B play a key role in many applications. In 
contrast to mg, they are determined by the QCD dynamics itself. We cannot yet 
calculate theoretically their values from first principles of QCD since it would require 
more or less exact solution of QCD in the strong coupling regime. Instead, we can 
try to measure them extracting from known properties of hadrons. 
The value of \j? g is known: since 

2mn 2 M 
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describes the interaction of the heavy quark spin with the light cloud and causes 
the hyperfine splitting between B and B*, 

/4 ~ - 2m b (M B * - M B ) ~ - (M%. - M§) « 0.36 GeV 2 . (71) 

In actual QCD y? G logarithmically depends on the normalization point; usual one- 
loop diagrams yield 

In the mass relation given above the operator is normalized at the scale (i ~ nib- 



J-Jtf. = -C G (,)^^ , C G c (^) ) 1 . (73) 
m Q 2m Q y a s (/i) / 

Evolving perturbatively to the normalization scale // ~ 1 GeV slightly enhances the 
value of but this effect is numerically insignificant, and we usually neglect it. 

The kinetic expectation value /^ 2 has not been directly measured yet. It enters 
various distributions in semileptonic decays or processes of the type b — > s + 7. 
First attempts to extract it from the semileptonic distributions were reported |28|, 
29| , however the analysis is vulnerable at some points, and so far the outcome is 
inconclusive. 

A model-independent lower bound was established @, 0, |], |J] 

fi 2 n > /4^0.4GeV 2 (74) 

which constrained possible values of /x 2 . One of its possible derivations will be 
given below in Sect. 4 (see also Sect. 5.1). Here we instead mention its physical 
interpretation. 

Unlike quantum-mechanical pure potential problems, in QCD the heavy quark 
kinetic operator is expressed in terms of the covariant derivatives ttj = iDj = idj + 
Aj which include the vector potential Aj. The latter lead to non-commutativity 
of different spatial components of the momentum operator in the presence of the 
chromomagnetic field B, 

2 

[vr^TTfc] = iG jk = -ie jM Bi ~ ——itjki ji ' 0-4 GeV 2 . (75) 

This non-commutativity immediately leads to the lower bound on the expectation 
value of 7? 2 [BUf, in full analogy with the uncertainty principle in quantum mechanics. 



Formally, it follows from the positivity of the Pauli Hamiltonian [31 



More physically, it means the Landau precession of a colored, i.e. "charged" particle 
in the (chromo)magnetic field, Fig. 3. Hence, one has (p 2 ) > \B\. Literally speaking, 
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in the B* meson the quantum-mechanical expectation value of the chromomagnetic 
field is suppressed, (B z ) = —[Iq/3. It completely vanishes in the B meson. However, 
the essentially non-classical nature of the 'commutator' B proportional to the spin 
operator of the light clound (e.g. (B 2 ) > 3(B) 2 ), in turn, enhances the bound which 
then takes the same form as in the external classical field. 




Figure 3: Landau precession of a free charged particle in the magnetic field. The average 
of the momentum square is bounded from below even in the absence of binding potential. 
This illustrates the physical meaning of the inequality between /j, 2 and [Iq. 

It is worth emphasizing that the inequality between the kinetic and chromomag- 
netic expectation values takes place for /z 2 normalized at any point /z, provided ji 2 G 
is normalized at the same point. For large /i it becomes uninformative; so, it is in 
our best interests to use it at fi = several units xAqcd ^ 1 GeV. 

The B meson average of the kinetic operator was estimated using the technique of 
the QCD sum rule in |32|, |33| and in |3^| , with quite different results (0.5±0.15) GeV 2 



and 0.1 ± 0.05 GeV 2 , respectively. A dedicated discussion of the subtleties inherent 
to those analyses can be found in Ref. [0], Sect. 6. The discreapancy roots to a 



different treatment of the contributions of the excited states which are a background 
in applications of the QCD sum rules, and the actual result seemingly lies somewhere 
in between the two extremes, probably closer to the upper value. It was later argued 



in Ref. [[35fl based on the analysis of the QCD sume rule approach to the three-point 
correlators in the harmonic oscillator, that using the same continuum threshold for 
the two- and three-point correlators (which is routinely assumed in QCD sum rules) 
leads to a systematic bias, and that it is advantageous to use a lower continuum 
threshold for the three-point functions. Applying this prescription to modify the 
analysis of Ball et ai, the authors got a value close to 0.4 GeV 2 . 

It is important to note, however that the normalization point dependence of the 
kinetic expectation value was not treated consistently. In fact, the above papers 
attempted to determine the HQET parameter — Ai rather than the properly defined 
QCD expectation value /i 2 dependent on the normalization point. In simple QM 
models with fixed numbers of constituents in a hadron, one would have them equal, 
— Ai = ji 2 . However, in actual QCD — Ai can be thought of as the kinetic expectation 
value /i 2 (/i) (/i is the normalization point of the operator) from which "all pertur- 
bative" (/i-dependent) pieces are subtracted. In this respect — Ai is a close relative 
of the pole mass, and likewise cannot be defined consistently in reality. Techni- 
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cally, it corresponds to extrapolating the normalization point in /i 2 (/-0 down to zero 
momentum, which leads to the usual conceptual problems. Related questions are 
discussed in more detail later in Sect. 4.2 and in Ref. @], Sect. 3. In practical terms, 
since the above mentioned QCD sum rules analyses were performed at the level of 
one-loop 0(a s ) perturbative corrections, it is more or less sufficient to use identifi- 
cation — Ai — > — Ai + ^f/U 2 — a4 (/■*)■ This suggests adding about 0.15 GeV 2 to the 
quoted values of — Aj to arrive at the QCD kinetic expectation value normalized at 
the canonical scale around 0.7 to 1 GeV. 

To summarize, the value of the kinetic expectation value /i^ in B mesons deter- 
mined from the QCD sum rules is about 

nl ~ 0.5 ± 0.15 GeV 2 , (76) 

however the literally quoted intrinsic uncertainty of the method may be underesti- 
mated here. It is noteworthy that this estimate seems to be in a good agreement 
with the bound Eq. (|74|). More details regarding various subtleties in the kinetic 
operator can be found in Ref. 0. 

3.3 Heavy quark potential 

We discuss here briefly some aspects related to the heavy quark potential in QCD. 
Strictly speaking, it refers to a somewhat different situation where two heavy quarks 
(actually, a quark and an antiquark) are present simultaneously. Heavy quark po- 
tential attracted recently attention in connection to pair production of ti and bb near 
the threshold, and in particular in view of its apparent connection to the problem 
of the heavy quark mass illustrated in Sect. 3.1. 

A closer look at the notion of the heavy quark potential reveals certain subleties. 
As a matter of fact, it is not completely clear what exactly must be called the 
potential between heavy quarks in actual QCD. 

The original notion of the potential refers to the interaction of infinitely heavy 
(static), or completely nonrelativistic heavy particles, which is instantaneous. The 
most familiar example is the electromagnetic interaction of heavy charges. The 
Hamiltonian of such a potential system is given by 

where m« are masses of particles and the potential V is a function of their instant 
coordinates. Taking the limit — > oo (at fixed r*j, which corresponds to semi- 
classically high excitation numbers of the quantum system in Eq. (|77j)) eliminates 
quantum uncertainties in the coordinates and allows to measure the potential di- 
rectly as the position-dependent energy of the infinitely slowly moving collection of 
particles. This is well known for QED where the potential of the charges q\ and q2 
(in units of electron charge) is given by 

^ ED (,R) = a em (0) ^ . (78) 
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This expression is exact in the absense of light charged particles; the known quantum 
corrections appear only if other matter fields are not much heavier than the scale 
1/R. 

The definition of the similar quantity in a non-Abelian theory like QCD is more 
tricky. First, one has to limit consideration to the systems in the colorless state. 
The heavy quark potential is defined only between Q and Q in the color singlet 
case; even there the requirement of gauge invariance is not trivial. The color of the 
individual heavy quark remains fully quantum in nature and changes through the 
interaction with gluons. In contrast to usual coordinates, the limit niq — > oo does 
not make color variables describing the state of the heavy quark semiclassical. 

To avoid this problem the heavy quark potential is defined via the vacuum ex- 
pectation value of the long Wilson loops: 



1 

T-^oo T X V Jc(R,T) 



V(R) = - Jim - In (Tr Vexp[if A^dx^ ) , (79) 



where the rectangular contour C(R, T) spans distance R and T in the space and time 
directions, respectively. (It is usually assumed to be in the Euclidean space as in 
Eq. (f79"P.) This definition is intuitively clear, since Wilson lines describe propagation 
of the infinitely heavy quarks. Moreover, such Wilson loops are readily computed in 
the U (1) gauge theory (free QED) and, of course, reproduce Eq. (|78|) (with qi = — q 2 ). 
In QCD the perturbative expansion of V(R) has been computed through order 
Usually the potential in the momentum representation is considered: 



VTO - / d* V {S) e -| ^ { l + (f - f „,) I + c 3 (^ + . . .) . 

(80) 

There are important peculiarities in thus defined interaction in the QQ system; 
they were first analyzed by Appelquist et al. already in the late 70's f37j. Due to 



gluon self-interaction, the potential in higher orders contains diagrams of the type 
shown in Figs. 4. Here the dashed line denotes the Coulomb quanta (they mediate 
instantaneous interaction in the physical Coulomb gauge), while the wavy lines are 
used for actual transverse gluons. These diagrams signal the real propagation in 
time of transverse gluons. This means that at this level the problem ceases to be a 
two-body one, and includes more full-fledged dynamical degrees of freedom. 

Moreover, while the first diagram Fig. 4a with the rung gluon appearing in order 
safely converges at the gluon momenta ~ 1/R, including additional Coulomb 
exchanges, as in Fig. 4b leads to infrared divergence. With one exchange it is 
logarithmic; the formal degree of the infrared divergence increases with adding extra 
Coulomb quanta between the emission and absorption of the transverse gluon. The 
physics behind this infrared behavior was discussed in Ref. ||37||: emission of the soft 



transverse gluon changes the overall color of the QQ pair and, therefore modifies 
the interaction energy between them. The energy shift associated with the exchange 
of the transverse gluon depends nonanalytically on the energy denominator (in the 
language of non-covariant time-ordered perturbation theory), since the gluon can 
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be arbitrarily soft. Formally expanding the exact result in a s includes expanding 
in this change of the Coulomb energy proportional to a s /R, therefore one obtains 
increasing infrared singularities. These arguments suggested that resummation of 
the Coulomb exchanges in Fig. 4b would render these diagrams finite, however the 
effective infrared cutoff is of the order of a s /R and, thus the potential is not infrared 
finite perturbatively containing terms ~ 1/R- hxa s starting order O (af). 
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Figure 4: Diagrams for heavy quark potential in QCD. Dashed lines are instantaneous 
Coulomb exchanges, transverse gluons propagating in time are shown by wavy lines. 

a) The convergent diagrams with the transverse gluon as a rung. 

b) Adding more Coulomb exchanges inside the ladder with the transverse gluon leads to 
infrared divergence in perturbation theory. 

This infrared singularity emerges since the external momentum scale 1/R does 
not fix all (spacelike) momenta in the problem. In the static situation the low- 
momentum tail in the propagation of very low-momentum (ultrasoft) gluons is cut 
off only by the interaction energy of the quarks, which appears itself only in the 
perturbative expansion; therefore it turns out to be suppressed only by powers of 
a s . As a result, beyond order the heavy quark potential, whatever it means, is 
not a short-distance quantity even at arbitrary small space separation between the 
heavy quarks. [] Of course, this becomes possible only since time interval used to 
define V(R) is taken infinitely large. Should one use instead finite-T Wilson loops, 
the time interval would place the final infrared cutoff in all gluon exchanges. 

This purely perturbative analysis shows that there is no direct analogue of the 
potential between heavy quarks in QCD. The QQ system incorporating all gluon 
interaction is not a two-body system but includes actual propagation of gluons with 
energy small compared to 1 /R. The interaction then cannot be universally described 
by an instantaneous potential and is intrinsically non-local in time. An attempt to 
integrate out in one way or another the "extra" degress of freedom does not yield 
in general a meaningful analogue of the potential: the result would not be universal 
but rather pro cess- dependent. Thus, the prescription (f[|) based on Wilson loops in 
QCD does not yield the heavy quark potential in its conventional understanding. 

5 This situation is not peculiar to non-Abelian gauge theories. Similar long-range propagation 
of light charged particles occurs also in QED if their masses are well below 1/R. 
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In practice, this may not pose serious problems in the perturbative computations. 
Through order all diagrams including Fig. 4a are convergent and saturated by 
the gluon momenta of order 1/R. The above nonlocalities in time are then revealed 
only at time intervals smaller than R. In many applications with nonrelativistic 
heavy quarks the scale 1/R is still much smaller than the characteristic dynamic 
energies. In the spirit of the Wilsonian OPE such exchanges then can be integrated 
out to yield a local in time effective potential. Yet it must be remembered that this 
can be done only up to a certain order in a s , depending on the concrete problem. 

As discussed above, the heavy quark potential V (from now one we assume that 
it is defined via Wilson loops, Eq. Q79D) is not an infrared finite object already in the 
perturbative expansion. The situation can be better for the integral of V(q) over 
spacelike momentum: 

v ^=Swf v(i!) - (81) 

In order to get a finite result we need to assume an ultraviolet cutoff at a certain scale 
/i, for example, adding the step-function 6(n 2 —q 2 ) to the integrand. It was mentioned 
in Sect. 3.1 that in electrodynamics such an integral of the potential equals minus 
double of the contribution to the mass of the static source coming from the same 
range of gluon momenta. In QCD the perturbative diagrams for both ttiq and V(q) 
are more complicated. Yet the same relation holds in order 0(a 2 ) as well. Most 
simply this is seen in the Coulomb gauge where all the effect to this order reduces 
to dressing the propagator of the Coulomb quanta. (An alternative discussion can 
be found in Ref. J38|.) There are reasons to believe that such a relation may hold 



to all orders in perturbation theory: the infrared contribution to V(0) in Eq. (plD 
equals minus twice the same contribution to the pole mass of a static source. 

Indeed, let us imagine we were able to introduce in some way an ansemble of 
gauge field configurations where the modes with momenta much larger than a certain 
scale \i are practically absent. This field-theoretic system would not need regular- 
ization, and everything can be expressed in terms of the bare parameters, including 
the bare quark mass mg'. At R— >0 the Wilson loop will approach the free value N c 
thus yielding V(0) = 0. This is clear on the physical grounds: QQ form a dipole with 
the infinitesimal dipole moment, and its interaction with any soft gluon field van- 
ishes as R goes to 0. It is important at this point that our gauge ansemble explicitly 
includes only soft modes. Otherwise, as in full QCD, the modes with \k\ ~ 1/R 
generate growing attractive potential at arbitrary small R. 

Next we note that V(R) is traditionally determined up to a constant; in pertur- 
bation theory (in four dimensions) the potential is actually defined as 

U(R) = V(R) - V(oo) ; (82) 

we assign U to this "standard" potential to distingush it from V(R) which has 
a precise meaning in a finite theory. While U(R) by definition vanishes at R — > 
oo, V(oo) does not and reflects nontrivial interaction with the gluon field. In the 
momentum representation V(q) contains explicitly the term V(oo) S 3 (q), which is 
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discarded in the usual perturbative computations. Since at R— > oo the Q and Q lines 
are well separated, their interaction (at least in perturbation theory) must vanish as 
1/R, and the value of the Wilson loop is simply given by the mass renormalization 
of each static source: 

V(oo) = 2 (mg cn - mg } ) = 26m Q . (83) 
Therefore, for the "ordinary" potential we have 

U(0) = / tSs K eg (?) = -2Sm Q . (84) 

Here V reg is the usually computed regular part of V(q) not containing self-energy 
diagrams yielding S 3 (q). 

The above relation for the infrared contributions to the mass and the potential 



0<?) = -25 m m Q (85) 



at first may seem strange: What happens to the perturbatively infrared-singular 

4 _»2 

contributions of Fig. 4b which, at order a*, behave like |f ln 2 ^, with e an infrared 
cutoff in the "rung" gluon momentum? The corresponding contributions are absent 
from the mass since the Coulomb exchanges are instantaneous. 

The subtlety resides in the necessity to introduce an upper citoff ji on gluon 
momenta in V lR . This not only cuts off V(q ) at q 2 ^> fi 2 , but also modifies V(q) at 
\q | ~ // and must lead to an additional term residing around \q | ~ \l which contains 

4 

a similar term %lne with the opposite sign. They do not, of course, eliminate 
dependence on e in the resulting ^ IR ((f), but have to cancel it when integrated over 
all q. 

This requirement has a transparent meaning: the amplitudes of emitting a very 
soft gluon out of the QQ system where initially Q and Q are set to be at arbitrarily 
small separation R, must vanish unless virtual gluons with k ~ 1/R are allowed to 
resolve Q and Q. Therefore, with all gluon fields softer than a certain /i, the sum of 
all diagrams for the emission amplitude (given by the part of the graphs in Figs. 4 
to the left from the transverse gluon) must vanish upon integrating over q . 

Relying on relation (^), one can try to define perturbatively a certain running 
heavy quark mass which is free from the leading renormalon uncertainty ~ Aqcd-Q 
One defines 



m P Q S ^ = < c +i / ^- 3 V Teg (q) = HdR V(R) ^ " R 

L J\q\<n {27l)° 71 JO 

(86) 

where simply the literal V(q) computed to a certain order in perturbation theory 
without an explicit cutoff is used, and the pole mass is taken to the same order in a s . 



R 



-\x cos fiR 



3 This idea was first discussed myself in 1996 and later independently advocated by M. Beneke 
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The mass mQ S (/j) is known as the "potential-subtracted" mass [58]. Although purely 
soft gluon configurations would generate 'an infrared' part of the potential V(q) 
literally different from V(q )9(\q \— p) used in Eq. (j86l) , this ansatz is hoped to include 
all (perturbative) infrared contributions originating from the domain sufficiently 
below /i, and simply subtracts some additional ad hoc perturbative pieces. 

A technical problem would arise here starting order - since V(q) becomes IR 
singular here, this routine literally fails. According to the previous discussion, the 
weight itself with which the infrared contributions are to be subtracted apparently 
gets modified by the explicit cutoff, should we implement it in practice. Therefore, 
they are not properly subtracted from the pole mass by prescription fl56|). It is con- 
ceivable that this mismatch exists already at order a^, but is not revealed explicitly 
due to infrared convergence of V to this order. 

It should be noted that V(0), although infrared finite in perturbation theory, 
is not a genuinely short-distance quantity (nor its /i-dependence) for which, for in- 
stance, usual OPE can be applied. Eq. (|86| ) also shows that at arbitrary large p it 
includes contribution from the potential V(R) at large distances R. The nonper- 
turbative definition of such an object therefore remains unclear. These interesting 
questions fall outside the scope of the present review. 



4 Heavy Quark Sum Rules and Exact Inequalities 
in the Static Limit mg^oo 

Discussing the basic parameters of the heavy quark theory and the properties of 
the actual b —>■ c transition it is often advantageous to resort to a theoretical limit 
when both quark masses are asymptotically large. This usually greatly simplifies the 
problem since eliminates proliferating effects of various power corrections in l/mq. 
In this limit physics of transitions between heavy quarks without change of the heavy 
hadron velocity is trivial: spin degrees of freedom are decoupled, only quasi-elastic 
transitions (those which do not change the light cloud, and the heavy quark spin 
can be flipped only by the weak current) occur between the respective members of 
the heavy quark symmetry multiplets. The transition amplitudes of heavy quark 
currents between the ground and excited states vanish in the heavy quark limit. 

Nontrivial physics emerges when the external weak current changes the velocity 
of the heavy flavor hadron. The most simple and instructive are transitions which 
occur in the lowest, linear in Av order. This so-called small velocity (SV) limit 
introduced by Shifman and Voloshin in the mid 80's as a theoretical tool for studying 
heavy flavor amplitudes, is also relevant for actual semileptonic decays of b to charm 
||39|| : the typical velocity of charmed final state hadrons is rather small. 

Physics of semileptonic transitions in the SV regime is described by the corre- 
sponding SV structure functions of the heavy mesons (or baryons; here we explicitly 



34 



formulate our consideration assuming the initial hadron is the ground state meson): 
W^(q ,q) = l.ha^-(H Q \ J d 3 xdx e^°*° iT{J,(x), Jj(0)}|ffg) (87) 

Q 

and we retain only terms through order v 2 . is normally expanded over invariant 
tensor structures. Heavy quark and SV limit lead to decrease in the number of 
structures for arbitrary weak currents - there is only two independent structure 
functions.^] We do not present here the corresponding formalism, but instead resort 
to an alternative consideration phrased in terms of individual transition amplitudes 
from the ground to the excited states. The analogies of the moments of the SV 
structure functions in this language are the so-called heavy quark sum rules. 

4.1 Sum rules 

The heavy quark sum rules we discuss are the following: 
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a sequence which, in principle, can be continued further. Here e& is the excitation 
energy of the fc-th intermediate state ( "P-wave states" in the quark- model language), 

e k = M H (k) - M Pq , 

7 Strictly speaking, the second one describing the antisymmetric Lorentz structure in the infinite 
mass limit would require considering the nonforward scattering amplitude with the overall change 
of velocity, and Hq carrying spin like the ground-state vector mesons. Alternatively, one can obtain 
its analogue considering the usual structure functions, for example, with vector currents to order 
1/toq even at zero recoil, or taking more complicated weak vertices with covariant derivatives. 
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while the functions r 3 ^ and describe the transition amplitudes of the ground 



state B meson to these intermediate states. We follow the notations of Ref. |f£C 

r(i/2) 



(F^|A^Q> = -Ti/a (vi-v 2 )^ , (96) 

and 



^-(^ 3/2) |^I^q) = -^iTs/ac^ (97) 

where 1/2 and 3/2 mark the quantum numbers of the light cloud in the intermediate 
states, j n = l/2 + and 3/2 + , respectively, and is the axial current. Furthermore, 
the slope parameter p 2 of the Isgur-Wise function is defined as 

' (P Q (v)\QloQ\P Q ) = l-Q 2 ^ + °(^) ■ (98) 



2M Pq ' vv ' v/ 2 



Equation fl88f ) is known as the Bjorken sum rule |4l| . Super convergent sum rules 



(|S9|) and (^T|) are new. Equation ([X]) was obtained by Voloshin [52]. The expression 
for /i 2 is the BGSUV sum rule [Q. The next one was derived in Ref. as well as 
Eq. (|95|). The last two sum rules are obtained along the same lines. The sum rule 
for the Darwin term p 3 D was first presented in . 



Let us illustrate the simple derivation of the sum rules for /i 2 and \Xq. Since 
in the heavy quark limit QQ pairs are not produced, we will use the quantum- 
mechanical language with respect to Q (but not the light cloud, of course). Moreover, 
as discussed in Sect. 1.1 it is convenient, at the first stage to assume Q to be spinless. 
The Q spin effects are trivially included later. Then the lowest-lying states, the S- 
wave configurations corresponding to B and B*, are spin-1/2 fermions, with two spin 
orientations of the light cloud. We shall denote them |f2 ); the spinor wavefunction 
of this state is \l/o- It is obvious that 

2 2 

(Oo|0(^;)(*A)Q|^o> = ^SjMVo - ^VlajiVo = J2(^\n) (k\7n\n ) , (99) 

6 b k 

where a complete set of intermediate states is inserted. They are spin-1/2 states (of 
the opposite parity with respect to \£lo)), generically denoted by <f>( n \ and spin-3/2 
states x ■ We will use the Rarita-Schwinger wavefunctions for the latter, i.e. a 
set of three spinors xi obeying the constraint <7iXi = 0- The normalization of these 
spinors is fixed by the sum over polarizations A 

£x<(A)x](A) = 4i-^ • (100) 

A 6 

Defining the reduced matrix elements a n and b m as 

(0("V^o> = a^a^o , (x (m) |7r;l^o) = & m xf %o , (101) 
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where (j)^ and stand for the states as well as for their wavefunctions, we get 

/4 = -6^|a n | 2 + 2^|6 m | 2 , (102) 



(71 



and 

£ = 3£|a n | 2 + 2]T|& m | 2 . (103) 



ni 



These expressions can be immediately generalized to the actual case of the spin-| 
quarks Q. The quantities a n and b m are to be understood as the matrix elements of 
biDb between the B meson and higher even-parity states. They are related to r| 

and t^m as follows: 



/2 



(n) a n (m) 1 f~tr\A\ 

T l/2 = — , r 3/2 = -ft— , ( 104 ) 



and, therefore, 



A = 3 fe^l 2 + E4|r$| 2 ) , ^ = 3 f2Ed|r 3 ( / m 2 ) | 2 - 2£e 2 Jr$| 2 ) . 

\m n / \ m n / 

Note that in atomic physics the combinations analogous to e |r| 2 are called "oscillator 
strengths" . 



Relations ( 104 ) are most easily obtained from the fact that the SV amplitudes 
of transitions to the corresponding states are given by the overlap (n(v)\B(v = 0)) 
(in the case of spinless Q; for spin-| they are given by (n(v )\Bq\B{v — 0))). Then 
we use the general relation 

\H Q {v)) = | J ff Q (0))+7r -V7f| J ff Q (0)) + 0(v 2 ). (105) 

which nicely elucidates the meaning of the small velocity sum rules: the operator 
n^^vTi) acting on \Hq) is the generator of the boost along direction v. Indeed, to get 
\Hq(v)) one must find the eigenstate of the Hamiltonian with heavy quark moving 
with the momentum q = rriQV. The only part which explicitly depends on momentum 

-*2 

comes from the heavy quark Hamiltonian (plus, in general, the higher terms 
in l/mq). We use the relation exp (— iqx ) Hq exp(zgx) = H,Q + v7r + mQV 2 /2 and 
drop the last term which is a constant; A obviously commutes with x. Then 
Eq. ( |105| ) represents the first-order perturbation theory in 57i = vn; the unperturbed 
Hamiltonian is Ho = ii (further details can be found in [[|, Eq. (178) and Sect. VI). 
In the second-quantized notations the very same relation takes the form 

\H Q (v)) = \H Q (0)) + J d i xQ^ 1 v^Q{x)\H Q {Q)) + 0(v 2 ). (106) 

The first four sum rules are obtained differently, using the OPE for the heavy 
quark transition operator. The Bjorken sum rule basically states that the probability 
of the quark to hadronize into some heavy flavor state after weak interaction, is 



exactly unity. In the different notations it was stated already in Ref. [jjq| . The 
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'optical' sum rule for A follows from the fact that no independent dimension-4 heavy 
quark operator exists for the forward matrix elements. This means that the apparent 
change in the kinematics when passing from the quark level to actual B mesons must 
be compensated in the average energy of the final states by the sums over excited 

-»2 

states: the threshold for quark transition is at go = ^— while for the hadron the 

-.2 -*2 

elastic peak is at go — 2 m h ■> l° wer by A^-- (The original derivation by Voloshin 
was different.) 

The new sum rules (|89D and ( PTJ ) can be obtained applying the OPE to the 
nonforward scattering amplitude where the momenta flowing in and out the two 
weak vertices are not equal, so that the overall change in velocity v = V\ — Vi is 
nonzero. They actually require the part antisymmetric in v\ and V2-f\ The new 
parameter of the heavy quark theory £ is unfortunately unknown. It is defined as 
a small velocity elastic transition matrix element between the states with explicit 
spin of light degrees of freedom: for the ground-state vector mesons like B* 

^^{HQ{v,e')\QiDjQ(0)\HQ(0, e)) = Vj(e'*e)+^^e'*v)ej—E'?(ev)}+0(v^ . 

(107) 

Alternatively, for spinless heavy quarks Q this would read as 

(fio(^)l WQ(0)|ft(0)> = ~ v 3 tfj^o - i| ew v k $1*1*0 + O (v 2 ) , (108) 

etc. Let us note that the full (quasielastic) matrix element in the l.h.s. of Eq. (|107|) is 
not well defined due to ultraviolet divergences in the static theory, and therefore the 
symmetric part proportional to A depends on regularization. Yet the antisymmetric 
part proportional to S is well defined and finite. It can be directly measured, for 
example, on the lattices. We expect £ to be about 0.25 GeV. 

Let us briefly outline the unified derivation of the sum rules Eqs. (|88p- fl9lD . The 
OPE approach leading to the sum rules is described in Sect. 5. Here, however, we 
consider the nonforward scattering amplitude similar to Eq. ( |S^ ) with = JJ = 
QloQ assuming that the final hadron has nonzero momentum Mh q v : 

T(q ;V,v) = ^±-(H Q (3)\ J d 3 xdx e'"*° iT{J (x), 4(0)}\H Q (0)) (109) 

with q = itiqv, and we retain only terms through second order in v and v. All 
terms suppressed by powers of 1/m.Q are discarded. As the energy variable we take 

e = go — ( \j <l 2 + m 'Q— m o) — go — mQ 2 ■ With this choice e = would correspond to 

8 This complication is not necessary for the sum rules for operators dimension 5 and higher 
(^G' PlS' etc.) where the antisymmetric structure functions appear, for example, in B* for the 
spacelike vector current vertices, although they are I/rag suppressed. However, the excited states 
contribute to these structure functions with extra factor of e 2 , and thus yield only the higher sum 
rules for local operators. 
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the elastic transition for free heavy quark; the larger e, the larger is the mass of the 
intermediate state. 

As explained in Sect. 5, we expand the amplitude Eq. ( |109| ) in 1/e assuming 
mQ^>eS>AQCD; the dispersion relation equates the terms in the 1/e expansion with 
the moments of the discontinuity of T due to on-shell transitions into intermediate 
states. The latter to the second order in velocities are the ground state (B^) and 
the P-wave excitations. The OPE yields 

-T(e;tr^) = i-^<flg(^)|g(l-?+?)Q(0)|ffg(0)> 
e ZM Hq 

+ \7^(H Q (v)\Q(tDv)Q(0)\H Q (0)) 
e zm Hq 

+ ^^^(H Q (d)\Q[( l Dv) 2 ^D ( l Dv)]Q(0)\H Q (0)) + ... (110) 



e 



Comparing this to the explicit hadronic contributions to the small velocity moments 
1/27T / dee k ImT(e; v, v) we get the stated sum rules. The structures antisymmetric 
in v and v — v emerge if Hq carries spin correlated with the spin of light cloud, 
and for B* yield the new sum rules involving the difference between the |- and 
^-contributions. The quasielastic transitions do not contribute to this structure for 
actual spin- 1 quarks. 

The computation proceeds similarly and is even simpler for spinless heavy quarks. 
The difference is only for the zeroth moment, sum rule fl39|): for scalar Q the anti- 
symmetric part is absent from the OPE expression for Jo, but the elastic contribution 
yields it with the opposite sign. 

The fact that the symmetric part of the D = 4 nonforward elastic matrix element 
(the first term in the r.h.s. of Eqs. ( |107| )-( |TU£| )) is given by ^ follows from equations 
of motion: 



M Hq (v-u)^(H q (v)\QQ(0)\H q (0)) = -(Wm) (H q (v)\QQ(x)\H q (0)) 

(H Q (v)\m Q QQ(0) - Q(z/VmM0)|#q(0)> 

(u^ = (1, 0, 0, 0) is the restframe four-velocity and Q are the original full QCD fields), 
or, to order v 2 

(M HQ -m Q ) ^ (H q \QQ(0)\Hq) = -v k {H Q (v)\QiD k Q(0)\H Q (0)) (111) 

which fixes the symmetric term in the above matrix elements. 

The sum rule Eq. ([59]) deserves a special note. The sum of the differences between 
t|/ 2 and Ty 2 amounts to a nonvanishing number which does not depend on strong 
dynamics. This can naturally raise a question of how this can be true, for the | and 
\ states are differentiated only by spin-orbital interaction. For example, the light 
quark in the meson can be, in principle, arbitrary nonrelativistic as well, in which 
limit this interaction can be switched off. 
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The point is that in the nonrelativistic case r 2 are large scaling like inverse square 
of the typical velocity of the light quark, r 2 ~ l/^ S p, and the relativistic spin-orbital 
effects must appear at the relative level ~ v 2 p simply because spin ceases to commute 
with momentum to this accuracy. This then leads to corrections of order 1 in the first 
two sum rules Eqs. (|88|) , (|89|) . In the Bjorken sum rule the relativistic constant 1/4 
is obscured by perturbative effects, while the sum rule Eq. (^) shows this universal 
relativistic effect in the most clean way 

As is seen from the preceding analysis, physics of the SV transitions is described 
by two "structure functions" W + and W_: 

W + (e,v) = 2^E|r 3 ( / m 2 ) r5(6-6 m ) + .; 2 Ek 1 ( / 1| 2 5( e - e n) + ^V( e + A|i), 

m n 

W.(e) = 2j2\T$\ 2 S(e-e m ) - 2 £ |r$| 2 <5( e -6 n ) (112) 

m n 

and the corresponding invariant hadronic tensors 

m c m c n c n c ti — 

| r ( m )|2 W^l 2 

h_(e) = 2 -2 ^^Z 2 !. (H3) 

They directly enter, for example, 1/m corrections to various weak currents at zero 
recoil. Sum rules Eqs. (|89|) and (|9ll) represent the exact low-energy theorems for 

fc_(e). 

Let us note that the sum rule (|89D provides the rationale for the fact that pseu- 
doscalar mesons B, D are lighter than their hyperfine partners B*, D*. If the sum 
rule for fj,Q is dominated by the low-lying states then [Aq must be of the same sign 
as the constant in Eq. QS^D, which dictates the sign of the hyperfine mass splitting. 

The sum rules (p8|)-(p4]) obviously entail a set of exact QCD inequalities. They 
are similar to those which have been with us since the early eighties ^6f , and reflect 



the most general features of QCD (like the vector-like nature of the quark-gluon 
interaction). The advent of the heavy quark theory paved the way to a totally new 
class of inequalities among the fundamental parameters. As with the old ones, they 
are based on the equations of motion of QCD and certain positivity properties. All 
technical details of the derivation are different, however, as well as the sphere of 
applications. The first in the series is the Bjorken inequality g 2 > 1/4 We, 
in fact, have a stronger bound p 2 > 3/4. Other bounds include p 2 > liq and 
Pd ^ |pisl/2, Pd > ~Pls an d A > 2S which are valid at arbitrary normalization 
point. We also have the direct inequalities between the parameters of different 
dimension: 

2 3A 2 3 3 A 3 3 (p 2 ) 3 / 2 , ^ 
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All these inequalities are saturated if only one excited state contributes. A number 
of additional inequalities involving the energy of the first P-wave excitation Ai is 
discussed in Ref. |[21|| . For further reference, let us also write an obvious consequence 



of the two equations for /4 and /i^.: 

t n I ' 1/2 



/4-/4 = 9E4l^| 2 . (us) 



4.2 Hard QCD and normalization point dependence 

The sum rules (|88|)-(|94"|) express heavy quark parameters, including A = M B —rrib, /4 
and fiQ as the sum of observable quantities, products of the hadron mass differences 
and transition probabilities. The observable quantities are scale-independent. How 
then, say, A = M# — m?,, /i^ and [jl 2 g happen to be //-dependent? 

The answer is that in actual quantum field theory like QCD the sums over excited 
states are generally ultraviolet divergent when efc ^> Aqcd; in contrast to ordinary 
quantum mechanics they are not saturated by a few lowest states with contributions 
fading out fast in magnitude with the excitation number. The contributions of 
hadronic states with e^. ^> Aqcd are dual to what we calculate in perturbation 
theory using its basic objects, quarks and gluons. The latter yield the continuous 
spectrum and can be evaluated perturbatively using isolated quasifree heavy quarks 
as the initial state. The final states are heavy quarks and a certain number of gluons 
and light quarks. It is the difference between the actual hadronic and quark-gluon 
transitions that resides at low excitation energies. 

Therefore, in order to make the sum rules meaningful, we must cut off the sums at 
some energy // which then makes the expectation values //-dependent. The simplest 
way is merely to extend the sum only up to < /i, and this is the convention 
we normally use. Thus, in reality all the sums in the relations (|88| ) -(|94l ) must 
include the condition e& < //, which we omitted there for the sake of simplicity, and 
all the heavy quark parameters are normalized at the scale /t. The exception are 
super convergent spin-nonsinglet sum rules fl8§|) and (PT|) where in the perturbative 
domain /a <C Aqcd such //-dependence is power suppressed by factors c^qcd//* 
and can be neglected. For analytic computations it is often convenient to apply 
the exponential cutoff factor e _<Efc ^, which is essentially the Borel transform of the 
related correlation functions. Since at large // the cutoff factors differ only in the 
perturbative domain, the difference between various renormalization schemes can 
be calculated perturbatively. 

The high-energy tail of the transitions to order a s is given by the quark diagrams 
in Figs. 5 with 

in n J ^ 

where (u, k) is the momentum of the real gluon. The spin-singlet amplitudes are 
just a constant proportional to g s , and performing the simple calculations we arrive 
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at the first-order term in the evolution of fJ%([i)' 

d/4(M) 4 a 



d/t 2 3 7i 



+ ... . (116) 



Purely perturbatively, the continuum analogies of T1/2 and r 3 / 2 are equal and a 
similar additive renormalization of /ig and p| 5 is absent. 



\ 1 





Figure 5: Perturbative diagrams determining the high-energy asymptotics of the heavy 
quark transition amplitudes and renormalization of the local operators. 



The perturbatively obtained evolution equations ( |116| ), 
mine the asymptotic values of T1/2 and T3/2 at e ^> Aqcd: 



sE^I^ + E^' 2 



a, e 



1/2 



E 



e ml T 3/2 



(m) 1 2 



-E4k ( 



2 |^(n) 1 2 
/2l 



9 7T 

3a s (e) de 
e 



ede 

-^E 



allow one to deter- 



(117) 



2 1 Ml 2 \p 2| («)|2l niS^I 
2^ A I ^ -m\ T 3/2 I C nl r l/2l ' ' 1 1 1 ;> ; 



ifm<£ 



£n<£ 



The last bracket is simply |//£,(e). Eq. (|117| ) can be extended to higher orders in a s , 



this simply reduces to using the dipole coupling a^(e) discussed in Sect. 2.4.1 

V _ 13 W 

2 4 J 71 



^ ( e - 5 /3+ln2^ 



+ 0(a 



(119) 



(a s is the standard MS strong coupling). Therefore we get a number of exact 
perturbative evolution equations [IP 



dp 2 (/i) 




d/i 


9 71 


dA» 


16 aW(/i) 


dp 


9 7T 
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(120) 
(121) 
(122) 



(the evolution for the Darwin operator is already more complicated since includes 
the usual anomalous dimension |4]]] ) . The similar questions for the alternative com- 
bination of the oscillator strengths entering /z^ and p 3 LS has not been studied yet. 



42 



4.3 On the saturation of the sum rules 



The question of the saturation of the heavy quark sum rules (in particular, the 
lower ones Eqs. fl88D - fl93"D ) is of primary importance for phenomenology of the heavy 
quark expansion. The sum rules state what is the asymptotic value in the r.h.s. 
at large enough normalization point /x, and perturbation theory tells us their /j 
dependence. From which scale /Jo this behavior applies, is a dynamical question. 
For super convergent sum rules ( |89| ) and fl9~T|) this is the question at which scale the 
sums approach the stated values with a reasonable accuracy. In order to sensibly 
apply quantitative 1/rriQ expansion, one must have rriQ > /i , possibly, rriQ ^> /i . 
While this is probably the case for b particles, such a hierarchy is not obvious a 
priori in charm. 

The existing numerical evaluations of A and fi 2 at the scale near 1 GeV sug- 
gest rather large values, approximately 0.7 GeV and 0.6 GeV 2 , respectively, which 
impose rather tight constraints. These facts are often neglected under various pre- 
texts, including challenging the accuracy of the QCD sum rules evaluations of /i 2 . 
The significant uncertainties, in principle, cannot be excluded. In a certain respect 
determination of Mb — rrib from e + e~ — > bb is somewhat indirect as well including 
mild theoretical assumptions. Indeed, even though QCD ensures that the quark 
mass appearing in the analysis of B decays and bb annihilation is the same, extrap- 
olating from the decay (scattering) kinematics appropriate for B decays with small 
negative q 2 for the b^b transitions, to the timelike domain q 2 ~ 4m 2 in the anni- 
hilation channel involves the long path. It brings in significant perturbative effects 
which are to be accounted for. If actual dynamics of gluonic degrees of freedom sets 
in to the perturbative regime late in energy, the associated uncertainties may, in 
principle, increase. 

Nevertheless, we point out that certain constraints following from the sum rules 
can be confined within the world with a single heavy quark, and they still are 
tight. Namely, the value of n 2 G ~ 0.4 GeV 2 has barely been challenged as extracted 
almost directly from B^ and D^*' masses. By virtue of the sum rules, the value 
of /i 2 is at least as large. Thus, regardless of the accuracy in evaluations of the 
kinetic expectation value, the question can be phrased in terms of the generally 
accepted value of /i^. At which minimal scale /io the value of /^(/io) reaches 0.3 
or 0.4 GeV 2 ? If this scale is below 1 GeV, large A and /j, 2 are almost inevitable. If, 
however, [Iq(1 GeV) is significantly below 0.4 GeV 2 , the chances for success in I/tuq 
expansion in charm are slim. 

The estimates for r 3 / 2 and Tin for the lowest P wave states group around 0.4, with 

4/2 ^ 4/2 - 400 to 500 MeV ( for the 

review see ^S|). Similar values were reportedly 
extracted from the overall experimental yield of the corresponding charmed P wave 
states f4"9||. It is evident that such oscillator strengths fall short in saturating the 



sum rules: 



6$ 2 /4 ~ 0.2 GeV 2 , 4/2 A - °- 3 GeV 



In principle, they alone should not necessarily, even though the idea of the dominance 
of the lowest states contribution is very appealing. Let us mention that in the 
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't Hooft model all the heavy quark sum rules are saturated with amazing accuracy 



by the first excitations [50, [51]. Is such a possibility excluded in QCD? Probably, not 
completely. The dominance of the first excitation with e ~ 500 MeV (recall that one 
must use the asymptotic uiq — > oo values of the excitation energies and amplitudes) 
is still possible if QCD sum rules underestimate the value of r^-Q Experimental 
determinations of r's are also questionable since l/m c corrections are not accounted 
for there. The estimates in the 't Hooft model suggest that they can be very large. 
In the cases when they are known explicitly in QCD, the l/m c terms generally turn 
out to be very significant as well [[52] . 



Another - and, apparently, the most natural - option is that there are new states 
with the masses around 700 MeV with similar, or even larger Ty 2 ~ 0.4 to 0.5; they 
can be broad and not identified with clear-cut resonances. The ^-states must be yet 
depleted up to this scale. All such states can be produced in semileptonic b decays 
and observed as populating the domain of hadronic invariant mass below or around 
3 GeV. It will be important to explore these questions in experiment. 

The similar constraints and the pattern of saturation follow from the sum rule 
(pl|) and, in particular, (p9[). The most natural "favorable" solution would imply 
existence of | states around e ~ 700 MeV with significant r 3 / 2 ~ 0.5, and suppressed 
transition amplitudes for the \ states in the whole domain of e below 1 GeV. 



5 Heavy Flavor Sum Rules; Finite rriQ 

In practical applications we usually need to know the decay amplitudes for actual 
b particles; therefore, the mass of the initial heavy quark, although large compared 
to Aqcd, is finite and the corresponding corrections must be taken into account. 
The c quark in the final state is only marginally heavy, and these corrections are 
typically very significant. In b — > u transitions the final state quark is light and 
the literal 1 /mq expansion becomes inapplicable. However, in certain cases - for 
example, in the sum rules - the energy of the final hadronic state rather than m q 
itself is an expansion parameter, and one still can use the similar expansion as long 
as the transferred spacelike momentum is large Q. 

The finite-mg sum rules are less compact for a number of reasons. First, one 
needs to account for the various 1/mg-suppressed terms; yet the corrections are ex- 
pressed in terms of the hadronic parameters we have encountered already analyzing 
the static SV sum rules. This is true if one expresses everything in terms of the 
quark-level kinematics. Since in practice we count all the energies from the observ- 
able hadronic thresholds, the nonlocal correlators which determine the expansion of 
the heavy hadron mass in powers of I/tuq enter as well. This is the only way how, 
for example, A enters the sum rules. Another source of nonlocal correlators of heavy 

9 Let us note that the technology of the QCD sum rules assumes the approximate duality starting 
e = 1 GeV or even lower (the energy are counted from the heavy quark mass there) . Therefore 
accepting poor saturation of the exact heavy quark sum rules at this scale and still relying on the 
QCD sum rules predictions is not self-consistent. 
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quark operators is the deviation of the expectation values over the initial hadron (B 
meson) from their asymptotic values which would exist if mf, — > oo. 

Second, beyond the static approximation external weak currents proliferate, and 
as much as five independent structure functions describe the decays even for the 
standard V— A current. Likewise, we often have to go beyond the SV approximation 
and do not assume that the final state velocity is small. Therefore, the structure 
functions as in usual deep inelastic scattering (DIS) become dependent on two kine- 
matic variables, say go an d q 2 ■ However, to keep the parallel with the heavy quark 
and the SV limit it is often advantageous to choose instead go and |g | = \Jq 2 — q 2 - 
(This choice is also convenient in a number of practical applications to b — > c decays.) 

The sum rules are derived in QCD using the standard methods of the short- 
distance expansion [|J . One starts the analysis from the forward transition amplitude 

T (12 \q -q) = — L J d 3 x dx e l " 0X0 (B\iT{cT w b(x), bT {2) c(0)}\B) (123) 

where I^ 1 ' 2 ) are some spin structures or, more generally, local operators. The tran- 
sition amplitude contains a lot of information about the decay probabilities. As 
usual in QCD, one cannot calculate it completely in the physical domain of q. The 
amplitude ( |123| ) has several cuts corresponding to different physical processes. The 
discontinuity at the physical cut g < M B — <J M^ + q 2 describes the inclusive decay 
probabilities at a given energy released into the final hadronic system. The cut 
continues further than the domain accessible in the actual decays, see Fig. 6. 
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Figure 6: Cuts of the transition amplitude in the complex go plane. The physical cut for 
the weak decay starts at go = Mg — (M^+g 2 ) 1 ' 2 and continues towards go = — oo. Other 
physical processes generate cuts starting near go = ±(M^ c +g' 2 ) 1 / 2 (one pair); another pair 
of cuts originates at close values of go- An additional channel opens at go>2mb+m c . 

For the practical case of left-handed currents = 7^(1 — 75), r^ 2 ) =7^(1 — 75) 
and the hadronic tensor T^ 12 )(go; q) can be decomposed into five covariants |53| 

= -hg^ + h 2 v il v v - ihse^apVnqp + h A q^q u + h 5 (q^v u + q u v^) , (124) 

each having the same analytic and unitarity properties. In particular, they obey the 
dispersion relation 

h Mo) = 77- / U hjkj!*k + polynomial , (125) 

27T J On-Oo 



45 



where ^ = (1,0,0,0) is the four- velocity of the decaying meson and Wi are observable 
structure functions: 

w l = 2lmh i . (126) 

For decays into massless leptons the structure functions W4 and W5 do not contribute 
due to conservation of the lepton current; they are proportional to mf and are 
relevant only for the decays into r lepton |54|| . 

In QCD we calculate the amplitude ( |123| ) away from all its cuts. Essentially, it 
can be expanded in the inverse powers of the distance from the physical cut, e : 



e = M B -^/M 2 D + q 2 -q = mt-^+g 2 + <K<f )-<?o , (127) 

where 

S(q 2 ) = (M B -m b ) - ^M 2 D +q 2 - ~ O ( A Qcd||t) + O . 

The "off-shellness" e must be chosen in such a way that that |e| and |e-arg e] ^> Aqcd 



but, simultaneously, |e| <C yml + q 2 , m b - The second requirement allows us to 
"resolve" the contributions of the separate cuts. 

How does the deep Euclidean expansion of T^ 12 ) help to constrain the amplitude 
in the physical domain of real e lying just on the cut? The dispersion relations 



Eq. ( |125| ) tells us that the coefficients in the expansion of the amplitude in powers 
of 1/e are given by the corresponding moments of the spectral density in T^ 12 \ 

= 1 rljaT^wm^, = _i £ ^ , ImT a V;?W; 

7r J e' — e 7r e fc+i J 

(128) 

the similar expansions hold for separate invariant structures. 

On the other hand, we can build the large-e expansion of the transition amplitude 
per se, treating it as the propagation of the virtual heavy quark submerged into a 
soft medium. The expansion takes the general form 



T {12 \q ;q) = 2^-/ d 3 fdx e^- 



(B\b(x)rV> 



{m c +ip) 



m 2 c -(iD) 2 -\<jG 



r( 2 ^(0)|S) 



x'O 



(r«K+m&7o-^]^r^ w-^ 9 — ^-r^ r( 2 )) 

m z c — (m b —qo)+q— 2m b ii + 2qTC — 7c z — ^aG 



(r^[m c + (E c +e-5(q 2 )) l0 + (n+ fix 

y [2(E c + e-5(q 2 ))n + 2q^+n 2 + -aG i 
h [-(e-5(q 2 m2E c +e-5(q 2 ))f +1) 
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Here E c — ^Jm^ + q 2 , the symbol [...] x o denotes the Ox matrix element of the corre- 
sponding operator (it can be understood as the integral operator in the coordinate 
representation), and we use the short-hand notation 



(•••) = ^f B j ^xdx Q (B\b{x)[..Ub{0)\B) ■ 

In the first equation we used the full QCD fields and then passed to the low- 
energy fields according to Eqs. fll^p-flU]). Picking up the corresponding term l/e k+1 
in the expansion over 1/t at e > Aqcd and evaluating the expectation value of the 
resulting local operators (e.g., 6(0)7r M 6(a;) = 8 A (x) 6^6(0)), one gets the sum rules 
sought for. Taking k = yields the sum rule for the equal-time commutator of the 
currents cT^b and bT^c; k — 1 selects the commutator for the time derivative, etc. 

Applying this general procedure to the heavy quark transition amplitude, we face 
a subtlety which is specific to the OPE in heavy quarks. Namely, speaking formally 
the dispersion relation represents the amplitude as a Cauchy integral over all cuts it 
has, while we try to reconstruct the structure functions (or their moments) referring 
to a particular cut describing the heavy quark decay process. It was first discussed in 
Ref. Performing the OPE we identify the integrals over the separate cuts in the 
quark Green functions with the corresponding cuts in the physical amplitudes. This 
identification is established order by order in 1/ttiq expansion since the contributions 
of the physical cut behave as powers in 1/e while the contributions of the other, 
"distant" cuts behave like l/(e + M) n where M ~ mq is the distance to the other 
cut, and therefore they are regular at e^rriQ. This property was given the name 
of global duality [|J. Beyond the practical version of the OPE global duality is, in 
general, violated and may play the role of additional theoretical background together 
with local duality violations.^ 

The problem of separating the extra cut contributions would become serious, 
for example, in the b — > u transitions if y/cp approaches m&. The u-channel cut 
then becomes too close to the physical decay cut |55| , and their contributions in the 



correlator Eq. ( |123| ) cannot be taken apart. This complication originates from the 
presence of the reverse ordering of the two weak decay vertices in the correlator. 
Naively, it could have been avoided if one considered the retarded - instead of the 
time-ordered - product in the correlation function, or simply defined the hadronic 
tensor as the dispersion integral only over the discontinuity of the decay cut (con- 
tinued towards larger q Q accessible, in principle, in the process of scattering of the 
weak current on the heavy hadron). However, in order to apply the technique of the 
OPE, we must have two necessary ingredients: first, the amplitude must have the 
proper analytic properties, and, second, be represented as the correlator of local field 
operators for which we can write the functional representation in the background 
gluon field. The combination of the two requirements dictates using the T-ordered 
product in Eq. (|123|) , and thus makes the question of global duality unavoidable. 

10 It is probable, however, that global duality is not too sensitive to the peculiarities of 
Minkowskian kinematics, and its violations are truly exponentially suppressed at large rriQ. 
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We quote here the expressions for the first three moments of the structure func- 
tions ifi.2,3 for the standard V — A weak currents [ffl; the higher moments require 
accounting for the nonperturbative effects beyond C(Aq CD ). The variable e de- 
scribing the excess of the final state hadronic energy above the threshold is defined 

as 

e = M B - ^M 2 D * + q 2 - q , (130) 

i.e., is counted from the energy of the D* meson (rather from the lowest state D 
contributing to the sum rules); this convention does not affect the zeroth moments. 
Moments of the structure functions are 



Z7T J 



e Wi{€, q- 



(131) 



where i labels the structure function; they can be considered separately for axial 
current (AA), vector current (VV) and for the interference of the two (AV). 
For the zeroth moments we have 
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Expressions for the VV functions are obtained from the axial ones by replacing m c — > 
— m c ; the structure functions n/ 1 ' 2 )^ and w^ AA,vv vanish. The above equations are 
analogies of the Bjorken sum rule; however they incorporate nonperturbative effects 
which appear at I/ttlq level; the corrections are not universal and differ explicitly 
for different currents and structure functions. 
The first moments look as follows: 
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At q = these relations determine I/tuq terms in the masses of heavy mesons. 
Their derivatives with respect to q 2 near zero recoil give the Voloshin's "optical" 
sum rule. Here they are obtained with better accuracy for arbitrary, not necessary 
small, velocity and incorporate 1/toq relative corrections. The latter appear to be 
quite sizable when q is not particularly large. 

The third sum rules, which are relations for the second moments of the structure 
functions, are calculated only in the leading non-trivial approximation. They look 
rather simple and manifestly satisfy the heavy quark symmetry relation 



2E C T (1)AA _ %Ec T (1)VV _ E c (2)AA _ E c (2)VV _ ow T (3)AV _ 
1=T~ 1 2 — 1 2 — 1 2 — 1 2 — —<m>c 1 2 — 

E c +m c E c —m c m;, mj, 

.£S^ + *^y. (138) 

All higher moments vanish in our approximation. We used above the notation Ejj* 
for the energy of D* and E c for the energy of the c quark in the free quark decay: 



E D * = ^M^+q 2 , E c = ^ml + q 2 . (139) 

The quantity (MB-^nb)—{ED*—E c ) which enters first and second moments determines 
the difference in kinematics between the free quark and the quasielastic B — > D* 
transitions. At zero recoil q = it is of order Aq CD , however in the general situation 
scales like Aq CD : 

(M B -m b ) - {E D ,-E C ) ~ 



The explicit expressions for the structure functions themselves including these 
nonperturbative corrections, and the semileptonic decay distributions are given in 
Ref . ]53| . Pure perturbative corrections to the heavy quark structure functions were 



computed in Ref. . 



5.1 Zero recoil sum rules; \V c b\ from B — ► D*lv at zero recoil. 
5.1.1 Sum rules for fi 2 and y? G 

Considering the semileptonic transitions driven by the pseudoscalar weak current 
J 5 = J d 3 xci'j 5 b(x) (i.e., at zero recoil q = 0) one obtains the sum rule 

1 r"...(5),^ J, _ V- IP 12 _ ( 1 



27T 



jy^e) * = £ \W = (^-^) 04M-/4M) (141) 



yielding the inequality fJ^(fi) > this is the field-theoretic analogue of the 

relation Eq. (|115|) . This correspondence is transparent keeping in mind the nonrel- 
ativistic expansion of the pseudoscalar current cij 5 b ~ (l/2m c — l/2m b ) ip^aTiipb at 
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zero momentum transfer. Only the transitions to the P-wave states survive here in 
the leading in I/tuq approximation. 

The sum rule ( p.02| ) for n 2 G can also be readily obtained in this way. For example, 
one can consider the antisymmetric (with respect to i and j) part of the correlator 
of the vector currents (r^ =7», r^ 2 ^ =7j); such a sum rule must be considered for 
B* mesons (in B the expectation value of P c h r vanishes). Alternatively, the sum 
rule for the correlator of the nonrelativistic currents ciTjb and b&irc directly can be 
considered. The OPE guarantees that all such relations are equivalent. 



5.1.2 F D * at zero recoil 

The concept of the heavy quark symmetry was very important for the evolution of 
studies of heavy flavor hadrons. The fact of the fixed normalization of the B — > 
D and B — > D* formfactors at zero recoil in the limit — > oo was of special 
significance in applications since suggested the method for determination of \V^\ 
from the B — > D* semileptonic decay channel near zero recoil. To this end one 
measures its differential rate, extrapolates to the point of zero recoil and gets the 
quantity \ V c b Fp* (0)|, where Fp* is the axial B — > D* iv formfactor. Since the charm 
quark is only marginally heavy, it is very important to estimate the corrections, in 
particular, nonperturbative. The exclusive transition amplitudes are not genuinely 
short-distance, such transitions proceed in time intervals ~1/Aqcd- Nevertheless, it 
turns out that the large-distance effects appear in these kinematics only suppressed 
by l/m 2 c b : 

F D *(0) = i + o(j£)+ Sf /m2 + 5f /m3 + ... (142) 

The absence of 1/ itlq corrections was first noted in Ref. B -H It can be readily 
understood. Let us consider, for example, the vector B^D transition at zero recoil: 

(D\cry b\B) = 2Jm b M d {\ + — - — + ...) (143) 

(short- distance effects are neglected). The relative magnitude of l/m.& and l/m c 
terms is fixed since at m& = m c all corrections must vanish identically. T-invariance, 
however, says that the coefficients for l/m c and terms must be equal since B 
differs from D by only the value of the heavy quark mass. Thus, both terms must 
vanish. Actually, this statement is the heavy-quark analogue of the Ademollo-Gatto 
theorem for the SU(3)fi breaking effects which is routinely exploited in determina- 
tions of \V US \ from K — > tt£u and semileptonic hyperon decays. This observation 
was later studied in more detail in Ref. 0] and is usually called Luke's theorem. It 
improves the credibility of this method of determination of \V c b\ in spite of certain 
experimental difficulties. 

Since the corrections to the heavy quark limit are governed by the mass of 
the charm quark, even the I/w/q effects a priori can be significant. The need in 

11 The first preprint version of the paper |37f discussing the heavy quark symmetry in QCD was 
given to me by M. Shifman in July 1986. 
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evaluation of the 1/ttiq corrections in Eq. ( |142 ) for practical purposes was realized 



quite early [58] . In these days the theory of the power corrections in heavy quarks 
was immature, so that it was hard to decide even the sign of <5jy m 2- Quantitative 
application of the heavy quark symmetry to charm was generally viewed overly 
optimistic; it was believed that the deviations from the symmetry limit in Fp* must 
be very small, at the scale of 2%, and suggestions that they could be as large as 
10% were categorically refuted [Q. At present, with the application of dynamic 
methods in the heavy quark expansion, the actual estimates of Fjj* rather fall close 
to 0.9 §. 

The existing estimates of the power nonperturbative corrections in Fp* are based 
on the sum rules for heavy flavor transition amplitude discussed in the previous 
section. Studying the zero- recoil transition we fix the spacelike momentum q — 0. 
The axial current 07^756 produces D*, Dir and higher excitations in semileptonic B 
decays at zero recoil. A straightforward derivation yields the following sum rule for 
this current (cf. Eq. (Pg) for 3h^ A at q = 0): 



D* 



£ |F,| 2 = fcM - Af /m2 - A? /m3 +0 — , (144) 

0<ei<fj, \ m Q/ 



where 




Af /- 2 = 3^ + - I — + —2 + oZTZT ) (145) 



(the l/m,Q corrections to the sum rule A^ m3 are also known, Eq. (8.21) of Ref. |2l|]), 
and the elastic formfactor F D * is defined as (D*(q = 0)\ 07^756] B) = y/2M B M D *F D »e*. 
Fi denote the axial-current transition formfactors to excited charm states % with the 
mass Mj = M^.+ej, and I; a is a short-distance renormalization factor. Contributions 
from excitations with e higher than fi are dual to perturbative contributions and get 
lumped into the coefficient £a(aO of the unit operator, the first term in the right-hand 
side of Eq. 



The role of \i is thus two-fold: in the left-hand side it acts as an ultraviolet 
cutoff in the effective low-energy theory, and by the same token determines the 
normalization point for the local operators; simultaneously, it defines the infrared 
cutoff in the Wilson coefficients. 

The sum rule Eq. ( |144|) leads to the upper bound: 

|F D ,| 2 ^a(/i)-Af /m2 -A^ /m3 - £ |F,| 2 , 

0<e,:<At 



1 . , . ML_~_Ml 

]>rni 



6? /m2 > -Af /m2 > B * B ~ 0.035 . (146) 



The last relation is a model- independent lower bound for the 1/m 2 corrections to 
Fd* at zero recoil 0. 

To obtain an actual estimate of the nonperturbative corrections rather than a 
bound, we need to know something about the contribution of the excited states in the 



51 



sum rule Eq. (144). Unfortunately, no model-independent answer to this question 



exists at present. The best we can do is to assume that the sum over the excited 
states is a fraction \ of the local term given by //^ and [1q, 



El^l 2 = X*fa, (147) 

Ei<IX 



where on general grounds x~ 1- The contribution of the continuum Dir state can be 
calculated 0, however theoretically it is expected to constitute only a small fraction 
of the sum over resonant states. Trying to be optimistic, we rather arbitrarily limit 
X by unity on the upper side, that is, put x = 0.5±0.5 @; the larger is the smaller 
is Fd*. 

Arguments in favor of the assumption \ < 1 are rather soft. In perturbation 
theory x = \ holds in the first order (see Sect. VII), but this relation changes in 
higher orders. The fraction \ was recently computed analytically |51[ in the exactly 



solvable 't Hooft model (QCD in (1+1) dimensions in the limit of large number of 
colors) and turned out to be close to 0.55; in this model it is almost saturated by 
the first radial excitation with e ~ 700 MeV. 

The short-distance renormalization factor £a(/-0 is calculated perturbatively. To 
the first order it was computed in 0, all-order BLM resummation performed in 
Pef . The most technically involved part of genuine 0(ot 2 s ) corrections to £a(a0 was 
computed in |5^| , and the complete result given in the non-BLM corrections 
turn out to be small slightly decreasing £a(a*)j with the overall estimate £a(p) — 
(0.99) 2 at jj, around 0.7 GeV. These computations, however were accomplished in 
expansion in n/rriQ and the terms to order /i 2 /mg were retained, which is consistent 
if only I/ttiq correction to the formfactor Fd* are addressed. In view of the presence 
of terms suppressed by powers of l/m c , the higher-order terms are expected to 
be significant. This was suggested by the perturbative BLM resummation itself 
where, in particular, \jvn? c corrections in ^ 2 were found to be at least at the scale 
of 2%. To estimate such effects, it is possible to evaluate £a(p) completely as a 
function of fi/rriQ, in the BLM approximation which presumably yields the dominant 
contribution. The typical numerical outcome of such computations is shown in 

1 /2 

Fig. 7. It is thus reasonable to accept the value £ A (ji) ~ 0.96 ±0.015 at /i^0.8 GeV 
for the short- distance renormalization of the axial current .0 The uncertainty here in 
the short-distance renormalization can hardly be significantly reduced further since 
the inherent momentum scale is relatively low. The 0(a s ) corrections to the Wilson 
coefficient of the kinetic operator in the sum rule Eq. ( |145| ) was also computed to 
the next-to-leading order in the correction turns out to be insignificant. 



12 A similar number, with much smaller uncertainty is often cited in the literature for the so-called 
ijA factor introduced in HQET to denote 'purely perturbative' scale-independent renormalization 

1/2 

of the axial current at zero recoil. Technically t\a coincides with £/ (/i) at /i — > to any order 
in the perturbative expansion. Regardless of intrinsic deficiencies of such a notion, the numerical 
value of t]a has the irreducible infrared renormalon uncertainty of at least 3-5% which by itself is 
much larger than the overall uncertainty in t]a quoted in the literature, see, e.g. [[n). The numerical 

1/2 

close coincidence of the value quoted for ija and (0.8 GeV) is accidental. 
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Figure 7: Short-distance renormalization of the zero recoil b— >c axial current as 

a function of fj,. AQCD= 250MeV. 



The 1/rriQ term in the zero recoil sum rule Eq. ( |144| ) -Af,3 is estimated to 



constitute about —3%. Although it further suppresses the apparent value of Fd*, we 
discard it in our estimates. It is expected to be only a smaller fraction of the overall 
1 / uiq effects in the B—> D* formfactor, and keeping A^ m3 is hardly legitimate since 
the leading correction S^ m2 in the formfactor is not known completely. 
Assembling all pieces together we get for \ — 0.5 ± 0.5 

Fd* ^ 0.89 - 0.015 ± 0.025 cxcit ± 0.015 pcrt ± 0.025 1/m 3 . (148) 

Estimates of the uncertainties in the contributions from higher excitations and of 
the magnitude of the (9(l/mg) corrections are not very firm and are rather on the 
optimistic side; they can be larger. With I/tuq corrections amounting to ~ 8% one 
expects |#i/ m 3 | > 0.025 simply on the dimensional grounds. 

5.1.3 Quantum-mechanical interpretation 

We already illustrated the quantum-mechanical meaning of the inequality between 
fi 2 and \Xq in Sect. 3.2. Let us look at the sum rule Eq. ( |144| ) from the similar 
perspective [|J. From the point of view of light cloud in B meson the semileptonic 
decay of the b quark is an instantaneous replacement of b by c quark. In ordinary 
quantum mechanics the overall probability of the produced state to hadronize to 
some final state is exactly unity, which is the first, leading term in the r.h.s. of 
Eq. ( |144| ). Why then are there any nonperturbative corrections in the sum rule? The 
answer is that the 'normalization' of the weak current 07^756 is not exactly unity and 
depends, in particular, on the external gluon field. This appears as presence of local 
higher- dimension operators in the current. Indeed, expressing the QCD current in 
terms of the nonrelativistic fields used in QM one has 



(aiD) 2 <Jk o k [oiD) 2 (aiD)a k (aiD) 
im 2 . 8m 2 4m c mb 



c lklb b <-> y?+ <| o k - 1 — + — — j— \+0 [ — )}<?„ 




(149) 
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The weak current 0757^6, according to Eq. ( |149| ) converts the initial wavefunction 
\l> b into \I>: 

* B ^ i = o k * a - { <^>fc + "f^- - («*W« g > +„.[,„ (15 0) 
Then it is easy to calculate the normalization of \& : 

The additional terms are just the nonperturbative correction in the right-hand side 
of the sum rule. It is worth noting that the first two I/rug terms in the bracket 
in Eq. (|149|) are the result of the Foldy-Wouthuysen transformation Eq. (]17D . Ad- 
ditional discussion, including the pedagogical illustration of the role of the latter 
transformation in the 1/w,q expansion, can be found in Ref. @. 



6 OPE for Inclusive Weak Decays 

Inclusive widths of the heavy flavor hadrons are examples of the genuine short- 
distance processes. The decays proceed at the space-time intervals ~ 1/ my, (more 
precisely, inverse energy release), and the widths are affected by the soft strong 
dynamics to the minimal extent. The decay widths are, perhaps, the best studied 
example of applying OPE to quantify dynamical effects in heavy quarks. Heavy 
quark symmetry per se is of little help here, in particular in b — > u or nonleptonic 
transitions. 

The general consideration of inclusive weak decays runs parallel to the treatment 
of a(e + e~^ hadrons). One describes the decay rate into an inclusive final state / in 
terms of the imaginary part of a forward scattering operator (the so-called transition 
operator) evaluated to second order in the weak interactions [[39| : 

lmf(Q^f^Q) = Im J d 4 x i T {C w (x), 4(0)} (152) 

where T denotes the time ordered product and C w is the relevant weak Lagrangian 
at the normalization point higher or about rrtQ. The space-time separation x in 
Eq. ( |152j ) is fixed by the inverse energy release. If the latter is sufficiently large in 
the decay, one can express the nonlocal operator product in Eq. (|152[) as an infinite 
sum of local heavy quark operators Oi of increasing dimensions. The width for 
Hq — > f is then obtained by averaging Im T over the heavy flavor hadron Hq , 

{H Q \2lmf(Q^f^Q)\H Q ) G 2 F m 5 Q \V KM \ V JUf)f {Hq^Hq^ 
— = T{H Q ^f)= ^3 ^c, (/x) — 

(153) 

with Vkm denoting the appropriate combination of the CKM parameters. A few 
comments are in order to elucidate the content of Eq. (|153|) . 
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(i) The parameter \x in Eq. ( |153| ) is the normalization point, indicating that we 
explicitly evolved from tuq down to fi. The effects of momenta below /x are lumped 

into the matrix elements of the operators Oj. 

( f) 

(ii) The coefficients c\ (/i) are dimensioned, they contain powers of 1/toq that 
go up with the dimension of the operator 0{. Using the normalization introduced 
in Eq. ( |153| ), one obtains on dimensional grounds 

-(/)/ \ 1 iTT \n\rr \ ml A QCD a s^ dl ~ 

c > 2mr a {H ^ H ^ ~ ^ 

with di denoting the dimension of operator Oj. The contribution from the operators 
with the lowest dimension obviously dominates in the limit vtiq— >oc. 

(iii) It seems natural then that the expansion of total rates can be given in powers 
of 1/rriQ. The master formula ( |153| ) holds for a host of different integrated heavy- 
flavor decays: semileptonic, nonleptonic and radiative transitions, CKM-favored or 
suppressed, etc. For semileptonic and nonleptonic decays, treated through order 
1/rriQ, it takes the following form: 



Sf) (H Q \QQ\H Q ) M (/) 1 (H Q \Q-^GQ\H Q ) M 
' 3 w 2M Hr , 5 w m 2 n 2M 



(154) 



where Wilson coefficients are of order unity. 

We pause here to make a few explanatory remarks on this particular expression. 
First, the main statement of the OPE is that there is no correction of order 1/mq 
||18|| . This is particularly noteworthy because the hadron masses, which control the 
phase space, do contain such a correction: M Hq =mq (l+A/mQ + 0(l/rriQ)^ ; the 

parameter A, different for different hadrons, does not enter the width. The reason 
for the absence of the 1/tuq correction in the total widths is twofold: the corrections 
to the expectation value of the leading QCD operator i s on ty ~ ^V 771 ^ an d 

there is no independent QCD operator of dimension 4 for forward matrix elements. 
Since the coefficients functions are purely short- distance, infrared effects neither can 
penetrate into them. 

A physically more illuminating way to think of the absence of corrections of order 
1/mQ is to realize that the bound-state effects in the initial state (mass shifts, etc.) 
do generate corrections of order 1/rriQ to the total width - as does hadronization 
in the final state. Yet local color symmetry demands that they cancel against each 
other, as can explicitly be demonstrated in simple models. It is worth realizing 
that this is a peculiar feature of QCD interactions - other dynamical realizations of 
strong confining forces would, generally, destroy the exact cancellation. A detailed 
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pedagogical discussion of physics behind this cancellation can be found in Ref. 0, 
Sect. 3.1. 

Second, the leading nonperturbative corrections are ~ 0(l/mq), i.e. small in 
the total decay rates for beauty hadrons. The first calculation of the leading non- 



perturbative corrections in the decays of heavy flavors was done in fig , |62| , |63| , 

Third, the four-quark operators (QTq)(qTQ) with different Lorentz and color 
structure depend explicitly on the light-quark flavors denoted by q. They, therefore, 
generate differences in the weak transition rates for the different hadrons of a given 
heavy flavor.^ They describe the effects of Weak Annihilation (WA) and Pauli 
Interference in B mesons, and Pauli Interference and Weak Scattering (WS) in 



heavy baryons. Their effects were calculated already in mid-eighties [47 



A note should be made regarding the subtlety in understanding the above four- 



quark expectation values |k5[ |66fl . They must be taken in the effective low energy 
(nonrelativistic) theory with respect to the heavy quark Q. In particular, the heavy 
quark field Q(x) must contain only the heavy quark annihilation operator, while 
Q(x) only the creation of the heavy quark. This differs from the full QCD fields 
which contain both annihilation of the heavy quark and creation of the antiquark 
in the Q(x) field (and likewise for Q(x)). The full-QCD four quark expectation 
values over, say, B meson {B\bY \q ql \b\B) include, strictly speaking, the interme- 
diate states with the bb pair (not related to the gluon conversion into bb), and not 
only light flavor hadrons. Even though the energy gap between the two classes of 
states is large, about 2rriQ, the latter contributions to the expectation values are not 
necessarily short-distance and rather governed by strong coupling dynamics when 
their energy is close to 2rriQ. This additional nonperturbative component is absent 
from the expectation values in the effective low-energy theory in the sector with 
a single heavy quark. Simultaneously, the analysis shows that such contributions 
are absent from the inclusive widths as well. Therefore, strictly speaking one can- 
not "match" in the usual sense the full-QCD expectation values onto those in the 
effective theory. The renormalization scale evolution of the four-fermion operators 
in the effective theory is governed by the so-called "hybrid" anomalous dimensions 
introduced by Shifman and Voloshin in the mid 80s j|7[] . 

Fourth, the short-distance coefficients c\ (/z) in practice are calculated in pertur- 
bation theory. However it is quite conceivable that certain nonperturbative effects 
arise also in the short- distance regime. They are believed to be rather small in 
beauty decays 



Fifth, a new matrix element appearing in OPE, not discussed so far, is the scalar 



heavy quark density. Its nonrelativistic expansion originally established in Ref. (18 
using the QCD equations of motion for the heavy quark field, follows from the 
identity ([16]) and takes the form 

(Hq\QQ\Hq) = (Hq\Q^ Q\Hq) - W V 2 J J W +0{l/m%) . (155) 



13 Expanding {Hq\Qi(jGQ\Hq) /ttiq also yields contributions of order l/m^; those are, however, 



practically insensitive to the light quark flavors. 
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Since (Hq\Qj Q\Hq) = 2M Hq , the spectator ansatz indeed emerges as the asymp- 
totic scenario universal for all types of hadrons, and holds up to 1/ttjq corrections. In 
addition to Qtt 2 Q, the second dimension-five operator is the chromomagnetic oper- 
ator QioGQ. Since QD 2 Q is not a Lorentz scalar, it does not appear independently 
in Eq. (|154j) . For the pseudoscalar mesons, for example, we have 

1 2 3 M 2 M 2 

(p e ieQip g) = i - Ar + 3^V^ + ° (1/m « ) • (156) 



The reason for the kinetic operator term to appear is quite transparent. The 
first two quantities on the right-hand side of the equation represent the mean value 
of the factor \J\ — v 2 reflecting the time dilation which slows down the decay of the 



quark Q moving inside Hq |j4|, [65|, |68 



Equation ( |155| ) is readily obtained in the heavy quark expansion if one uses 
proper nonrelativistic heavy quark spinors incorporating the Foldy-Wouthuysen 
transformation; in the context of HQET this procedure was advocated in Ref. |p9 |, 



but was largely ignored for a few years. 

Finally, Eqs. (|154j) -( |155|) show that the two dimension-five operators do produce 



differences in B versus A^/S^ versus fib decays of order 1/mq. To a small extent 
they can also differentiate B and B s via the SU(3) breaking in their expectation 
values. Differences in the transition rates inside the meson family are generated 
at order 1/mg by dimension-six four-quark operators. They are usually estimated 
in the vacuum saturation approximation which - although cannot be exact - rep- 
resents a reasonable starting approximation. There is an intriguing way to check 
factorization experimentally fl65 |: similar four-fermion operators enter semileptonic 



b — > u transition rates. Moreover, in the heavy quark limit the four-fermion op- 
erators populate mainly the transitions into the hadronic states with low energy, 
and thus show up, for example, in the end-point domain of the lepton spectrum 
where their relative effect is enhanced. More precisely, these effects are present at 
large values of lepton invariant mass q 2 and thus can be isolated in the cleanest way 
studying double differential semileptonic distributions. Considering the difference 
of the decay characteristics of the charged and neutral B's in this domain, one can 
measure these matrix elements and even feel their scale dependence. Further details 
regarding the "flavor-dependent" preasymptotic effects can be found in dedicated 
papers 0, 0. 

It is important to keep in mind that the OPE approach discussed in this sec- 
tion implies that all decay channels induced by a given term in the short distance 



Lagrangian in Eq. (|1 52|) are included. It is not enough to consider the states that 
appear at the free quark level. The final state interactions can annihilate, for exam- 
ple, the cc quark pair into light hadrons; electromagnetic interaction, if considered, 
can do the same. Disregarding the channels that can emerge due to such final state 
interactions can violate the general theorems: even 1/tbq terms can appear in such 
incomplete "inclusive" widths. For example, the correction to the b — > s + 7 width 
does have nonperturbative corrections scaling like 1/mb due to the effect of weak 
annihilation (in mesons, or weak scattering in baryons) of the light quarks with 
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Figure 8: a) Spectator effect leading to l/mj correction in the decay width 6^5 + 7. 
b) A different cut of the same diagram leads to an electromagnetic correction in the 
hadronic decay width. Terms l/m& and 1/ttt| cancel out in the sum of two decay widths. 
The solid dot in the vertex denotes the penguin-induced b — > sqq interaction. 



emission of a hard photon in the penguin-induced weak decay b — > sqq, see Fig. 8a. 
This effect would only cancel against the (virtual) electromagnetic correction to the 
hadronic penguin- induced width, Fig. 8b. |72|, |65| . 

For the very same reason the spectator-independent nonperturbative effects in 
the actual b — > s + 7 width can also have in practice small corrections which scale 
like 1/mb [73]. The effective b — > s + 7 interaction is generated by loop diagrams 
and, including contributions from cc (««) pairs with momenta not large compared to 
mi, is not genuinely pointlike. This part of the decay interaction is properly treated 
considering the original b — > ccs (b — > uus) operators as the weak decay Lagrangian 
C w , and then truly inclusive widths for such an interaction have to include the 
probability of the usual decay b — > ccs. The latter is affected by electromagnetic 
interaction, again by virtue of the same diagram as for the b — > s + 7 decay, but with 
a different cut. Of course, in experiment - in contrast to the OPE - these processes 
are completely different and are taken separately. 

The OPE predictions for the nonperturbative effects in the inclusive width dif- 
ferences are rather nontrivial. Consequently, a number of questions are often raised. 
First, is this analysis operating in terms of quarks and gluons, applicable to actual 
b hadrons, with confinement leading to the drastic change of the physical spec- 
trum which consists only of the colorless hadrons but not quarks? Second, if this 
is correct, how to formally justify this OPE? Third, how to compute the leading 
nonperturbative corrections to the decay rates? And, finally, it is advantageous to 
have a transparent physical interpretation of the OPE machinery. 

The answer to the first question is definitely positive, although is not trivial; the 
validity of the OPE approach has been challenged more than once, even over the 
recent years. The raised criticism covered different aspects, but basically reduced to 
one main point: due to 'brutal confinement' there is no "duality" between the OPE 
computations operating in terms of quarks and gluons, and actual hadronic decay 
rates, either for nonleptonic, or for any type of the underlying quark transitions. We 
do not dwell here on discussing the reasons (usually poorly substantiated) behind 
such suggestions. This question was recently discussed in much detail in the context 
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of exactly solvable two-dimensional 't Hooft model where, in principle, all hadron 



masses and decay amplitudes can be computed. In Refs. [66|, |75| the analytic 
summation of the rates for open decay channels was performed in the expansion in 
powers of I/ttlq, and compared to the OPE predictions for the model. It was shown 
that to all orders in I/tuq the OPE series could be computed, they reproduced ex- 
actly the asymptotic expansion of the actual inclusive widths, for both semileptonic 
and nonleptonic decays. 

The physical interpretation of the main OPE result, the absence of the 1/rriQ 
corrections to the decay widths has been briefly mentioned above. A dedicated 
discussion can be found in Ref. 0, Sect. 3.1. In the following sections we briefly 
address the third and the second questions, respectively. 

6.1 Sample computation 

Here we illustrate the computation of the leading nonperturbative corrections in 
semileptonic, nonleptonic and b — > s + 7 -type decays. For simplicity, we do not go 
beyond order 1 /ml, and also neglect masses of the final state quarks and leptons; this 
will make expressions more compact. For semileptonic and nonleptonic decays we 
still will refer to one of the final state quarks as the charmed one. The corresponding 
decay Lagrangian is given (we omit the CKM factors) 

G 

£w = —j= 07^(1-75)6 -6^(1 -7 5 )w ; (157) 

this will describe the semileptonic decays as well if we switch off strong interaction 
of the du pair. We do not introduce here the color contraction schemes explicitly; 
color indices can be easily taken care of in the end. For b — > s + 7 we take the decay 
Lagrangian in the form 

A 
2 



£ w = -sia^F^(l+ l5 )b. (158) 





Figure 9: Diagrams describing the the transition amplitude in the external field for 
computing the b — > s + 7 (left) and nonleptonic or semileptonic widths (right). Inter- 
action with the gluon medium enters through Green functions for the final state quarks 
(intermediate quark lines), and through equations of motion for the initial b fields. 



The computation of the transition operator T in Eq. (|152|) is most simply done 
in the Fock-Schwinger (fixed point) gauge x^A^x) = (see Ref. |7£| for details). 
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The transition operator in the coordinate space is given by the products of the 
fermion and photon Green functions like b(x) G(x, 0)-D(x, 0) 6(0), or by the product 
of three fermion Green functions b(x) G(x, 0)-G(x, 0) -G(0, x) 6(0) for the two cases, 
respectively. For massless particles the Green functions are 



G(x, 0) 



iyt iXa.G a i3lpl5 



2vr 2 x 4 



]tt 2 x 2 



+ ... 



Dfj, u (x, 0) 



J fJ,U 



4ir 2 x 2 



(159) 



(a note of caution: in the fixed point gauge one should distinguish G(x, 0) from 
G(0, x) for quark Green functions), and the expansion of the gauge field is given by 



A^x) 



x u G ni {0) + ——x v x p (D p G Vfl ) + 



2 • 0! 



G above is the dual field strength, G 
Minkowski space. For b — > s + 7 we thus have 



[XV 



3-1! 



(160) 



G a/3 . Note that we work directly in 



T, 



MX 2 



2n 2 x e 



A 2 b(x) 



b(x)a i 



6 # 

7T 4 X 8 



+ 



2ir 2 x 4 



^TT 2 X 2 



^(1+75)6(0) 



2tt 4 x e 



G a /3l/3 + 



;i+7s) b(0) 



(161) 



where the summation over /i, v is straightforward. For the decays mediated by 
the four-fermion interaction we consider both semileptonic and nonleptonic cases 
simultaneously, simply introducing the factor £ which would switch off interaction 
of leptons with the gluon field, £ = 1 for b cud and £ = for b —>■ civ: 



^4-ferm = 4iG F K^H 



" 1 2tt 2 x 4 



SlX 2 X 2 



2vr 2 x 4 



^Tl 2 X 2 



7m 



4G 2 F b(x) 



k 2tt 2 x 4 

7T 6 X 10 



8tT 2 X 2 J 



7,(1-75)6(0) 



rGaplp + ... 



1-75) 6(0) 



(162) 



47T 6 X 8 

where we have used the Fiertz transformation for one of the V — A vertices, and 
the short-hand notation xGj denotes the repeating structure x a G a ^p. We observe 
that only the interaction with the antiquark is present to this order (using a variant 
of the Fiertz transformation, it is easy to show that this holds for arbitrary fermion 
masses [fL8| ). 

Now we need to make the Fourier transformation to pass to the momentum 
representation of the transition operator: 

r 2 



x 6 



71' 



— — p^ In {—p ) + polynomial in p 



7T 



— p 2 p,. In 1 
48 y 



-p 2 ) + polynomial in p 2 



7l~ 



X 



10 



64-24 



p 4 p Ai ln(— p 2 ) + polynomial in p 2 



(163) 
(164) 
(165) 
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(the divergent pieces come from infinite momenta running inside loops and can- 
not yield imaginary part). The polynomials do not have discontinuity and can be 
discarded for computing decay widths; it comes only from the logarithm. 

Now we need to evaluate the expectation values bp 4 $b, bp 2 p'b etc. which enter 
at p 2 To do this we recall that p M = iD^—A^ and in the Fock-Schwinger gauge 

one has 

^(0) = , D, t Au lo = ^GV(O) (166) 

(this implies, for example, that in our approximation [D^, A^] can be replaced by 
zero). Our strategy is to pull the gluon potential to the left since A(0) = and 
the result does not vanish only if a derivative acts on A^. In this way we arrive at 
the very simple rule that in this gauge 

bp 2n p'b = m 2 b n+1 bb - nmf' x b l -aGb + O (m^/mf) • (167) 

We also have an additional structure which, for the forward matrix elements reduces 
to the chromomagnetic operator using equations of motion for the b field: 

— ~ — % 

bp a G al3 -f^ 5 b = m b b -oG b . 



Recalling that T = 2ImT and collecting all terms, we arrive at [18, 33 




G 2 F m\ (B\bb\B) f /4 / 1 
sl 192tt 3 2M b \ m 2 \™l 

_ \ 2 ml (B\bb\B) f & 
Tb ^~ 2M~ B \ l - {l + l) ml + 

_ G 2 F m\N c (B\bb\B) \(c 2 + + c 2 _ c\-c 2 \ ( 2j%\ c\-c 2 _ 8^ ( 1 \\ 
nl 192tt3 2M B IV 2 2N C )[ m 2 J 2N C m 2+U \ml)y 

(168) 

Account for the charm quark mass is straightforward and yields the factor (1 — 
m c/ m b) 3 l° r the direct contribution of the operator /j,q, the second term in T n \. 
An additional piece expressed in terms of the derivative of the free quark phase 
space factor z with respect to m 2 comes from translating the result from the Fock- 
Schwinger gauge, cf. Eq. ( p. 671 ); it introduces the factor (1 — m^/m^) 4 for \j? G in the 
first term, and likewise in Y s \. The correction to the nonleptonic width b — > ccs with 
two heavy quarks in the final state was computed in Ref. |62|. 



6.2 How OPE can be justified for inclusive widths 

In the previous section we illustrated the calculation of the leading nonperturbative 
corrections to the inclusive decay widths. The expected question is: what such a 
computation of the quark decay width in the external gluon field has to do with 
the decay width of the actual heavy flavor hadron? The justification does not differ 
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conceptually from, say, the classical case of e + e annihilation to hadrons and makes 
use of the analytic properties of the transition amplitude and the dispersion relations. 



This was discussed in Refs. [65], |F?|, [74, ^6| with varying degree of detalization. 



To study the analytic properties of the forward transition amplitude we must 
introduce an auxiliary complex variable uj: 

A{uj) = J d*x (H Q \tT {£ w (x), £i(0)}\H Q ) (169) 

(once again, v is the four- velocity of the decaying heavy hadron). This a;-dependent 
amplitude can be visualized as the transition amplitude governing the total (weak) 
cross section of the scattering of a fictitious spurion particle 5* on the heavy quark, 

S(q) + H Q (p) light hadrons , (170) 

or the weak decay width in the process 

Q — > quarks (leptons) + S . (171) 

Such processes would appear if the weak decay Lagrangian is modified from, say the 
conventional four-fermion form to the "four-fermion + spurion" interaction, 

C w (x) -> S(x)C w (x) . (172) 



For simplicity it is convenient to assume, as in Eq. (|169|) that the spurion field does 
not carry spacelike momentum but only energy. 

The amplitude A(uj) has the usual analytic properties: it is analytic and has 
a number of cuts describing the physical processes in different channels (s, t or 
u). The physical cut corresponding to the weak decay of the heavy quark we are 
interested in, starts near u ~ E r where the energy release E r denotes tjiq minus the 
sum of the masses of the final state quarks and/or leptons. Other cuts are located 
far enough from this point and from the physical point uj = 0. The discontinuity 
across the physical cut at which the point uj = is located, describes the total decay 
width we are interested in. The OPE for the inclusive widths relies on the fact that 
the short- distance expansion of A{uj) runs in l/(u — E r ) and can be applied near 
the physical point uj = exactly as in e + e~ annihilation near a positive value of 
s^>Aq Cd . To the same extent, in principle, a certain smearing can be required if 
the hadronic probabilities still exhibit the resonance structure. 

Thus, there is no theoretical peculiarity in the asymptotic applications of the 
OPE, say for nonleptonic widths compared to semileptonic. It does not make a 
conceptual difference to perform a short- distance expansion of a single quark Green 
function (like in semileptonic widths or b — > s + 7), of the product of two Green 
functions (e + e~ annihilation) or of the product of three quark Green functions (the 
nonleptonic widths). 

Smearing in uj can often be phrased as smearing over the interval of rng. Indeed, 
in the heavy quark limit the decay amplitudes depend on just the combination my-oj, 
therefore 

A{uj, m Q ) ~ A(Q,mQ-uj) (173) 
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(there are power corrections to this relation associated with explicit mass effects in 
the initial state). Smearing in the heavy quark mass may look more transparent 
physically when the effect of opening new threshold is effectively averaged. 

The structure of the cuts in the amplitude A{uj) is particularly simple if all the 
masses of the final state quarks are large compared to Aqcd- What happens if some 
final state quarks are light, for example, even in semileptonic h —>■ u decays? The 
physical cut corresponding to the decay into u quark can come very close to the cut 
which would describe the scattering of u on the heavy quark - the distance between 
them is twice the energy of the quark, and it can be small if the heavy quark line 
is soft (i.e., both energy and spacelike momentum are small). We would not be 
able to distinguish between the contributions of the two processes if their respective 
cuts are too close in the scale of strong interactions. The answer is that in the 
OPE we remove explicitly all soft lines (including the quark ones) which carry small 
momenta; their contributions are described by the "condensates" where the soft legs 
are treated as the external fields acting on the initial state. The remaining pieces 
describe only energetic particles for which the proximity of the cuts cannot take 
place. From this perspective the OPE cutoff fi in computing the Wilson coefficients 
acts similar to assigning a mass fi to each potentially soft quark line. 

It is important to emphasize in this respect that the OPE treatment of these 
contributions in the "corners of the phase space" differs from, say, the computation of 
the regular 1/mg terms given by the chromomagnetic operator. In the latter case the 
effect (at least a part of it) comes from the kinematics where the corresponding final 
state quark is hard, and we really use the short- distance expansion of the propagator 
to compute this term. In the case of soft quark legs we do not compute them in the 
OPE, but, in a sense, parameterize their contribution. The large momentum flowing 
through the rest of the diagram only ensures that it is given by the expectation 
value of a local heavy quark operator. However, even if we knew the full transition 
amplitude with some accuracy near u = 0, we still would not be able to resolve 
the contribution of the physical and the "w-channel" cuts. These subtleties used to 
cause some confusion in the literature; the proximity of the cuts was viewed as the 
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6.3 \V c b\ from the total semileptonic B width 

Including the leading corrections, the semileptonic width has the following form 



ia 62, 33, 54 



^^\4-W- 2 ^i~l^ + ^ (174) 



where ellipses stand for higher order perturbative and/or power corrections, zq and 
di depend on m 2 c /ml. Regardless the exact value of the direct l/rn| corrections 
to r s i are rather small, about —5% and lead to the increase in the extracted value 
of \V c b\ by 2.5%. The higher order power corrections are at a percent level. There 
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is one new operator, the Darwin term which appears to order l/(m& — m c ) 3 ; it was 



evaluated in Refs. [77], |7S| and may decrease the width by up to 4%. The remaining 
l/mf corrections enter only through the l/m& terms in the expectation values of 
the kinetic and chromomagnetic operators in actual B mesons. They include the 
'spin-orbital' expectation value p\ s as a l/mj part of the full QCD operator b^aGb, 
and the nonlocal correlators introduced in Sect. 2.2. Their effect is minor and is 
included in the uncertainties of the chromomagnetic and kinetic expectation values. 

The leading-order 0(a s ) perturbative corrections are known from the QED cal- 
culations in muon decay ||79|| . The BLM corrections fl80 |, a part of the higher-order 



perturbative series associated only with running of a s in the first-order loop di- 
agrams, were calculated to all orders f8lf . Their impact appeared to be small if 



the width is expressed in terms of the properly defined short- distance masses |p2 
The remaining non-BLM ('genuine') two-loop a 2 s corrections were expected to be 
moderate if one uses the heavy quark masses normalized at the scale around 1 GeV 



They were evaluated in a series of papers by Czarnecki and Melnikov |p9| , |83 
and were confirmed to be small. At present, no significant theoretical uncertainty 
remains in the perturbative corrections to the semileptonic width. 

The running mass of the b quark is known with high precision, Sect. 3.1.2. To 
determine m c , we can rely on relation (P3|) for the difference between rrib and m c in 
terms of the hadron masses. It turns out that in this way the direct dependence of 
the semileptonic width on p z D almost cancels against the one coming indirectly from 
mfc — m c , therefore the width becomes sensitive only to the exact values of pi and 
P 3 - 



Evaluating the theoretical prediction we get 



1 1 0.105 ) V r B J I 0.1 GeV 2 / 



fl-0.01 ft "' ( :, G f ) Vl + 0.007^^V (175) 
V 50MeV J \ 0.1 GeV 3 / V ; 

The main theoretical uncertainty at the moment resides in the value of p 2 , which 
comes from constraining m^ — m c in the heavy quark expansion. It is worth noting 
that this is the only place where we relied on the expansion in l/m c . It is vulnerable 
to possible late onset of the 1/mg expansion which would show up as the large ex- 
pectation values of higher- dimension operators. Yet we note that the meson masses 
are expected to be the most robust observables; it is well known in ordinary quantum 
mechanics that eigenvalues of the Hamiltonian are more stable than wavefunctions 
and their overlaps. Nevertheless, it seems important to have an independent direct 
determination of m c to isolate this potential problem. 
Finally, we arrive at the model-independent evaluation 



V rh \ = 0.0412 



fBR(B 



\ 0.105 
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,(/4-0.5GeV 2 ^ 
V 2 

-3 



1 - 0-012 Q1Ge y 2 J • (1 ± 0.012 pcrt ± 0.01 mft ± 0.02) , (176) 



where the last error reflects rag and higher power corrections as well as possible 
deviations from local duality. 

6.4 r s i(6 — ► u) and determination of \V u b\ 

Similar to the treatment of T(B — > X c £v), it is straightforward to relate the value 
of \V u b\ to the total semileptonic width T(B — > X u iv). The dedicated analysis was 
performed in Ref. 



lTr , /BR(B°->X u &/)V/l-55ps\i 

|Vuf,| 0.00445 I o""002 J V t / ' ^'^ pert 

±0.03 mb ±0.015 nonpert ). 

(177) 

The dependence on /i 2 is practically absent here. The complete perturbative cor- 
rections are known analytically two two orders |j35|l , and again turn out to be small 



if one uses ra&(l GeV) as the input. The contribution at the level of a few percent 
can be expected from nonfactorizable expectation values of the four-quark opera- 
tors at order 1/mf. At this level it would be advantageous to measure the b — > u 
semileptonic decay width for the charged and neutral B mesons separately |l5 . 



The most direct way to disentangle b — ► ulv decays from hundred times more 
abundant b — > civ without tagging the secondary charm decay would be to study 
the invariant mass of hadrons in the final state, Mj: 

r(B^X£u), M] c =(j2P£ d r) ■ (178) 



For the free quark decay one has M\ ~ in b —>■ u and M x = m 2 c for the b — > c 
transitions. In actual decays the mass can take values exceeding m n and Mo, 
respectively. The increase in mass can originate both perturbatively if a hard gluon 
is emitted in the decay, or through soft bound-state or hadronization processes. 

The leading soft effects turn out to significantly modify the Mx distribution 
due to the effects of primordial 'motion' of the heavy quark caused by bound state 
dynamics in the initial b hadron. It was pointed out already in [64| that the QCD- 



based OPE automatically leads to an analogue of this "Fermi motion" which had 



been introduced phenomenologically long ago, first in p6| and then elaborated fur- 



ther to the status of a well-formulated model in [37]. A detailed description of the 



Fermi motion in the framework of the I/uiq expansion was later given in 



and in f89 ]; it goes outside the scope of the present review. Here we only briefly 
mention the basic facts. 

The "Fermi motion" in QCD has certain peculiarities which are absent in the 
phenomenological models [[3fj, |9(|. The distribution over the 'primordial' Fermi 
momentum F(p) is replaced by a certain distribution function F(x), where x < 1 
measures the momentum of the b quark in the units of Mfi-mj. It is similar to the 
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leading-twist distribution function of deep inelastic scattering (DIS) on usual light 
hadrons. First distinction is that F(x) is one-dimensional; one can define only the 
distribution over a certain projection of the momentum. 

Second, F(x) depends essentially on the final state quark mass (more exactly, on 
its velocity). While it is the same for b —>■ uiv and b — > s + 7 where the final quark 
is ultrarelativistic and is given by the light-cone distribution function like in DIS, 
it is rather different for c quark in b — * civ where it is closer to a nonrelativistic 
particle; in this case the light-cone distribution is replaced by the so-called temporal 
distribution function. 

Finally, the QCD F(x) is normalization-point dependent. This property is well 
known already from usual DIS; the evolution of the heavy quark distribution func- 
tion, where normalization point is well below itlq is, however, much stronger. 

The leading soft effects generating the Mx spectrum emerge due to same physics 
which gives rise to the Fermi motion. They are quite significant. For example, 
the average invariant mass square of hadrons in b — > u gets the nonperturbative 



correction ~ Aqcd ■ mb HH : 



(M 2 X ) 



nonpert 



7 

To 



m b (M B -m b ) +0 (A| CD ) ~ 1.5 GeV 2 . 



179) 



The perturbative corrections lead to (M x ) peTt ~ ^m 2 , however in B decays this 
increase is still smaller than through the nonperturbative effects. 




2.5 3 




2.5 3 



Figure 10: The integrated fraction of the b — > uiv events §>(Mx)- ( a ) : Dependence on 
mi, which is varied within ±50MeV. (b): Effect of variation of by ±0.2 GeV 2 . 

To quantify the QCD effects on the M x distribution let us introduce, following 
Refs. [0, ^] the fraction oib^u events with Mx below a certain cutoff mass M max : 

1 /-Afmax AY 

*< m -' = F(fW„ iMx m,- (180) 

$(0) = and $(Mg) = 1 hold regardless of any dynamics. The main question to 
theory is whether it can calculate accurately enough $(M max ) with M max <^ 1.6 GeV. 
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The dedicated analysis was carried out in , and the conclusion appeared to 

be quite optimistic: b — > ulv decays are not expected to populate the domain above 
M x = 1.6 GeV too significantly. The typical theoretical predictions for $(M max ) 
are shown in Fig. 10. 



The corresponding experimental analysis was first attempted by ALEPH |93 
and recently performed by DELPHI [|9*4"| . The Mx spectrum was found in a nice 
agreement with the above theoretical predictions, and the value of was 
found about 0.10, with the uncertainty of approximately 20 percentage points; the 
latter is still dominated by the experimental error bars. 

6.5 Summary on \V c b\ 

Here we give a brief summary of the two methods of extracting \V c b\- The most 
precise at the moment is the value obtained from T S \(B): Eqs. ( |175| ) and ( |176| ) with 
the central theoretical input values shown there lead to 

N,0M3.(l-C,,(»l). (181) 

The overall relative theoretical uncertainty in this result is 5 t h ^ 5%. With future 
refinements we can expect reducing the uncertainty down to 2%. 

The B — > D*lv zero-recoil rate also provides a good accuracy. The exact exper- 
imental status of the measurements extrapolated to q = is not completely clear at 
the moment. Until summer 2000 the results from the LEP experiments and from 
CLEO seemed to be in a good agreement, yielding F D * |V^| — 0.035 with a typical 
uncertainty of ±0.001 stat ± 0.002 syst . Using this reported average value and the lit- 
eral estimate Fd* — 0.89 only a bit lower value of \V c b\ — 0.0395 emerges. Adding 
only one (systematic) standard deviation to -Fd»|V^| would yield just the central 
value for \V^\ extracted from the total semileptonic width - even without varying 
the numbers within theoretical uncertainties. 

It is interesting that the central theoretical value in this method, according to 
Eq. (|148|) exhibits the dependence on /i^ similar in magnitude but opposite in sign 
to Eq. (Pip : 

0.0395. (l + 0.015 W -°^f } ) . (182) 

The theoretical uncertainty <5 t h here constitutes probably 5 t h ~ 6%, however this 
estimate relies on additional theoretical assumptions. It is not clear how it can be 
decreased. 

Recently CLEO came up with the new value for F D »\V ch \ = 0.0424 ± 0.0018 stat ± 
0.0019 syst |95|] which is significantly higher than both their old number and the 
current LEP average value 0.0345 ± .0007 sta t ± 0.0015 syst p5| , and now has larger 



error bars. Taking it literally, one even observes 'overshooting' - the central value for 
\V c b\ under the same theoretical assumptions would be even higher than the result 
from r s i(£>). Since the shift exceeds the previously quoted error bars, it may be 
premature to draw a definite conclusion at the moment. 
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In any case, it is remarkable that the values of \V c b\ that emerged from exploiting 
two theoretically complementary approaches are very close. The progress was not 
for free: it became possible only due to essential refinements of the theoretical 
tools in the last several years, which prompted us, in particular, that the zero- 
recoil B — ► D* formfactor Fd* is probably close to 0.9, significantly lower than 
previous expectations. The decrease in F D * and more accurate experimental data 
which became available shortly after, in summer 1994 reduced the gap between the 
exclusive and inclusive determinations of \V c b\. 

It must be noted, however, that most of the existing theoretical analyses essen- 
tially assume that the approximate duality between the actual hadronic amplitudes 
and the quark-gluon ones sets in already at the excitation energies ~ 0.7 to 1 GeV. 
In the formfactor analyses it is most obvious in applying the perturbative computa- 
tions to charm quarks and relying on l/m c expansion. In computation of the decay 
widths we, alternatively, invoke the local quark-hadron duality. While there are no 
experimental indications so far that this is not the case (at least, in the semileptonic 
physics), the proof is not known either. If that is not true, and duality generally 
starts only above 1 GeV, most probably one would have to abandon the idea of 
accurate determination of |V^,| from the exclusive B — ► transitions. The only 
option still open will be the inclusive semileptonic decays where the energy release is 
large, ~ 3.5 GeV. Of course, in such a pessimistic scenario the theoretical precision 
in \Vcb\ will hardly be better than 5%. 



7 Challenges in the HQE 

Concluding the review it is appropriate to mention cases where the heavy quark 
expansion apparently has problems. It is probably premature to speak of a direct 
contradiction to experiment; nevertheless, today's question marks carry the seeds 
of tomorrow's advances. Basically there are two problems where our theoretical 
understanding is lagging behind. Both are related to nonleptonic decays. 



7.1 Semileptonic branching ratio of B and T(b^ccs) 

The theoretical attitude to this problem changes with time. Twenty years ago the 
parton model gave a prediction BR s i(S) ~ 13-15%, which was accurate enough 
according to the existed standards. The variation reflected mostly the choice of 
quark masses. While the rates of c — > civ and b — > cud were affected by the choice 
of m c in more or less the same way, the rate of the b — > ccs channel decreased much 
faster when larger masses (closer to the 'constituent' rather than short-distance 
'current' ones) were used. Therefore, BR s i(I?) increases for larger masses. It should 
be remembered that the difference m^— m c is fixed, so the choices of m c and rrib are 
always correlated. Even though the heavy quark symmetry had not been formulated, 



the latter fact was clearly realized at least in the early 80's | 97| . 



The experimental situation was not very definite and indicated a rather large 
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BR s i(i?) which reasonably fitted the 'larger' mass option [ 9~Sfl . Since then it be- 
came standard to use the larger quark masses, and the parton model prediction 
for BR s i(£>) was accepted to be 13-14%, even though a smaller value could be ob- 
tained as well. Since the impact of the nonperturbative corrections was completely 
unknown and presumably significant, BR s i(i?) was not under intense scrutiny even 
when the better data became available. 

The situation changed when in 1992 Bigi et al. showed that there are no 1/rriQ 
nonperturbative corrections to the inclusive widths, both semileptonic and nonlep- 
tonic. The leading 1 jra\ nonperturbative effects were readily calculated |18|, j54] . 
They appeared to be suppressed, in particular, as a result of certain cancellations. 
The overall effect A npt BR sl (_B) was found to be about —0.5%; the nonperturbative 
effects could not be blamed for a discrepancy any more Hi. 



This prompted a more careful analysis of the perturbative corrections to the 
widths. In particular, the 0(a s ) corrections to the nonleptonic b — > cud width 
were calculated accounting for the nonzero m c ||100|| ; this additionally enhanced 



the nonleptonic width. Later the account for the charm mass was also completed 
for b — > ccs JTUTL p.02| , |103| |. Altogether, these 0(a s ) corrections further decreased 



BR s i(i?) down to 11-12%. Although this shift naively seems very significant and 
may raise concerns about the convergence of the perturbative corrections, it actually 
is not dramatic if one starts with more appropriate short- distance masses, the choice 
forgotten for historical rather than rational reasons. The relative increase in the 
nonleptonic width was particularly significant in the b — > ccs channel. 

The issue of the semileptonic branching ratio must be considered in conjunction 
with the charm yield n c , the number of charm states emerging from B decays. To 
measure n c one assigns charm multiplicity one to D, D s , A c and S c and two to 
charmonia. Zero is assigned to the charmless hadronic final state. It is obvious that 

n c ~ 1 + BR(E -> ccs q) - BR(5 -> no charm) . (183) 



Such a joint analysis was motivated already in p8[ : the energy release in b — ► ccs 
is rather small, and this can lead to significant duality-violating and higher-order 
effects. The stability of the perturbative expansion also downgrades. Measuring n c 
allows one to effectively exclude this channel from the theoretical calculations. 

At present the semileptonic fraction coming both from CLEO and LEP seems 
to be 1% to 1.5% lower than the 'preferred' theoretical expectation. The latter can 
be decreased playing with the ratio of the quark masses and boosting the effect of 
perturbative corrections - for example, taking strong coupling at a lower scale. This, 
however, increases n c beyond the experimental limits. In principle, allowing n c to 
lie in the upper corner of the experimental interval this would almost accommodate 
the LEP value, but still is somewhat off the CLEO interval for BR s i(I?) and n c . 

It is quite possible that this discrepancy should be taken seriously, and we would 
have to admit certain limitations in our ability to compute the nonleptonic B decay 
width. However, a more thorough analysis is badly needed which would be free 
from prejudices. The existing compilations rely on the old analysis where, among 
other things, pole masses of b and c are used as the input. This leads to a number of 
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problems: the perturbative corrections look too significant, and simply rewriting the 
first order result in different - but equivalent to order ct s - forms yields noticeably 
differing numbers. This, somewhat artificial, ambiguity is sometimes used to stretch 
the limits in a desired direction. Second, using the pole masses for quarks does not 
allow to immediately constrain their values which are nowadays well known in the 
proper scheme. 

Another place where refinement is needed, is a more accurate determination of 
the Wilson coefficient for /j,q in the semileptonic fraction. So far all the analyses 
used the expressions ( |168| ) from the original papers [18[ |62| where it was obtained in 



the leading logarithmic approximation only (barring the b — *> ccs channel, BR s i(5) 
is modified only via nonfactorizable diagrams appearing to order a s ), being propor- 
tional to a small product of two different color coefficients a\di of the weak nonlep- 
tonic Lagrangian. This may be a poor approximation, and the non-log piece coming 
from the extra gluon with k^rrib can dominate the result, as often happens with the 
LLA expressions in practice. Here there are special reasons to expect significance 
of such corrections: First, the corresponding diagrams can include gluon-gluon in- 
teraction, and such coupling is enhanced by color factors. Second, the LLA term 
suffers accidental cancellation between the different channels, the fact external to 
QCD itself. It is improbable that the chromomagnetic interaction itself can lower 
semileptonic fraction by, say, two percent. However, the existing apparent discrep- 
ancy between theory and experiment is lower now than it appeared eight years ago, 
and this certainly prompts a more careful analysis, to draw the final conclusion. 

7.2 Lifetimes of beauty hadrons 

As stated before, differences between meson and baryon decay widths arise already 
in order l/roq. The perturbative corrections to the lifetime ratios are completely 
absent. The lifetimes of the various mesons get differentiated effectively first in order 



l/m,Q. A detailed review can be found in ||104| , [70], [71 



Because the charm quark mass is not much larger than typical hadronic scale 
one can expect to make only semi-quantitative predictions on the charm lifetimes, 
in particular for the charm baryons. The agreement of the predictions with the data 
is good. I would even say it is too good keeping in mind that the l/m c expansion 
can hardly be justified. It is difficult to avoid mentioning one particular example. 
It is generally understood that the large lifetime ratio between charged and neutral 
D mesons, Td+/td° — 2.5 is due to destructive Pauli Interference in nonleptonic 
decays of D + . However, it is the semileptonic fraction of D + which is close to the 
'canonical' value of 20% obtained by simple-minded counting of the quark decay 
channels, while BR s i(D°) is much smaller. Yet, if one adds the l/m% shift in the 



semileptonic fraction due to the chromomagnetic term computed in Ref. [pC8| , the 
experimental values of the semileptonic fractions are reproduced. Clearly, one should 
not take such a coincidence too literally since the corrections of order unity discussed 
here cannot be rigorously justified in the framework of the standard 1/rriQ expansion 
itself. 
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Table 1: Predictions for Beauty Lifetimes (heavy quark expansion) 



Observable 


QCD Expectations (l/m& expansion) 


Ref. 


Data from |T05| 


r(B-)/r(B d ) 


l + 0.05(/ B /200MeV) 2 




72| 




1.070 ±0.02 


f(B s )/r(B d ) 


1± 0(0.01) 


1 


104 




0.945 ±0.039 


T(A b )/r(B d ) 


>0.9 




104 




0.794 ±0.053 



As far as the beauty lifetimes are concerned the 1 / m b expansion is to be applica- 
ble. Table [l] contains the recent data together with the predictions. The latter were 
actually made before data (or data of comparable sensitivity) became available. 

Data and predictions on the meson lifetimes are nontrivially consistent. Yet 
even so, a comment is in order for proper orientation. The numerical predictions 
were based on the assumption of factorization at a typical hadronic scale which is 
commonly taken as the one where a s (/ihadr) — 1- While there is no justification 
for factorization at fi ~ m&, there exists ample circumstantial evidence in favor of 
approximate factorization at a typical hadronic scale. More to the point, the validity 
of factorization can be probed in semileptonic decays of B mesons in an independent 
way, as was pointed out in |)5] . 



The possible effect of the nonfactorizable contribution has been discussed in 
detail in |65|, Q, and later in the dedicated paper Jfl]. They include the nonvalence 



gluon mechanism discussed long ago in [ 105 ] in the simplified language of the quark 



model. Significant nonfactorizable contributions would in general lead to large effects 
of Weak Annihilation in D s mesons where experimentally such effects are quite small. 
Of course, we cannot reliably treat in the I/tuq expansion, and Td s —T d + is 
sensitive to a particular combination of the expectation values of a few four-fermion 
operators. Nevertheless, this is an indication that the effects due to nonfactorizable 
expectation values must be suppressed. Later studies based on QCD sum rules and 
preliminary lattice results reproduced within intrinsic uncertainties the estimates 
one obtains relying on factorization at the low hadronic scale. The nonfactorizable 
effects may affect the difference between f(B s ) and r(B d ). 

The agreement of the data on B meson lifetimes with experiment is obscured by 
the apparent conflict for T\ b /rB d - To predict the l/rriq corrections to T\ b in the 1/rrib 
expansion one needs to evaluate the baryonic expectation values of two operators, 

(A b \bbu^ u\A b ) , (A b \bfbu-f ^u\A b ) . (184) 

They do not have the usual factorizable contribution, and their values are rather 
uncertain. Nevertheless, it was shown that their contributions cannot be too large 
[ 107f and the maximal effect in Y\ b does not exceed 10-12%. To achieve larger 



corrections one would have to go beyond the usual description of baryons where 
light quarks are "soft". The reason for such a limitation is rather simple. The 
four-fermion expectation values, in the language of ordinary quantum mechanics, 
are squares of the wavefunction at the zero separation between the heavy and the 
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corresponding light quark: 



*<to(0) = / T^^Mp) > ( 185 ) 

but, at the same time, the space integral of |^/Q g | 2 is limited by unity. A too large 
value of ^Q q (0) means that the wavefunction is sharply peaked at zero separation, 
which signifies presence of high-momentum modes. The above mentioned general 
bound agrees with the fact that the constituent quark model estimates typically 
yield only about 3 to 5% enhancement ||108|| . A similar conclusion has been reached 



by the authors of Ref. |[L09|| who analyzed the relevant baryonic matrix elements 
through QCD sum rules. A dedicated discussion of these questions can be found in 
Ref. 0. 



Recently, the preliminary lattice results for the relevant four-fermion expectation 
values were reported from the UKQCD collaboration, with the corresponding esti- 
mate ta 6 / T~B d — 0.90 lllOfl . We note here that the l/mn corrections to the lifetimes 



can well constitute 30 to 50% of the l/m,Q terms. Keeping in mind that the dom- 
inant effect in the A& lifetime is Weak Scattering, it is plausible that these effects 
further enhance the width. (This depends on the corresponding dimension-7 four- 
fermion operators with derivatives; taking the naive quark model picture one would 
simply replace in the Wilson coefficient by the square of the average diquark 
energy in A&, to account for these effects.) 

If future more accurate experimental data confirm the present deviation of T^ b /rB d 
from unity, and refinements of the theoretical estimates do not reveal unexpected 
enhancement, the most probable explanation of the discrepancy will be violation of 
local quark-hadron duality in nonleptonic decays. As a matter of fact, its significance 
in nonleptonic widths is theoretically expected a priori. Indeed, the expansion 
parameter for the widths is not mj, directly but rather the energy release which is 
noticeably smaller in b— >c. Moreover, the preasymptotic corrections depend on the 
concrete form of the weak interaction involved. For the four-fermion interaction they 
are enhanced: the arguments based on the high power of mass ttlq in the decay rate 
P^J suggest that the actual scale parameter is smaller, oc E re \/5. In the semileptonic 
decays this does not deteriorate the expansion since it is automatically protected by 
the heavy quark symmetry when m c increases. 

For the nonleptonic decays the heavy quark symmetry does not generally apply, 
and at insufficient energy release one expects significant violations of duality. The 
real problem here is that the few leading terms in the expansion have been evaluated 
and did not indicate the expansion blowing up at the level of 20%. 



8 Heavy Quark Expansion and Violations of Lo- 
cal Duality 

A large number of applications of the heavy quark theory are based on quark-hadron 
duality, more exactly, its local implementation. Although this notion becomes ex- 
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act at asymptotically high energies, at finite energy scale certain deviations must 
be present. How fast duality sets in and how large are these deviations are im- 
portant questions. These and similar questions are among most difficult. Yet they 
are important for phenomenology: determinations of \V c b\ we discussed rely on the 
assumption of local duality. 

The notion of duality in general terms was first introduced in the early days of 



QCD in Ref. [1 1 Ij but not pursued for quite some time. It was simply regarded as a 
problem of how the observables, say a(e + e~ — >hadrons) as a function of energy com- 
puted through quarks (and gluons) can be equated to the actual cross section. The 
former predict small smooth deviations from the would be free quarks production, 
while the latter has well-shaped resonant structure reflecting quarks permanently 
confined in colorless hadrons. Since such hadronization effects are nonperturbative 
in nature, the question of local duality is tightly related to treatment of strong in- 
teractions beyond perturbation theory. This was not really available in the early 
days of QCD; confinement was simply assumed not to affect averaged cross sections 
"significantly" . 

Nowadays we have in our disposal methods which are applicable - at least in 
principle - to quantify nonperturbative effects in a number of observables, including 
those of the actual Minkowski world. These are quantities amenable to study via 
the OPE. Local quark-hadron duality was given a new consideration a few years 
ago by M. Shifman |112| , who related its violations to the asymptotic nature of 
the power expansion in inverse large energy scale provided by the 'practical' OPE. 
These ideas were later reiterated and developed in a number of papers, in particular 



in connection to heavy quark physics (see, e.g., Refs. |74|j ). 



The basic observation can be illustrated on the example of heavy quark decay 
widths. The "practical" OPE yields the width in the power expansion 

££« = j4o + A + 4 + ... . (186) 

Tq m Q m z Q 

If the series in 1/mq were convergent (to the actual ratio), Th q would have been 
an analytic function of ttlq above a certain mass mo pointing to the onset of the 
exact local parton-hadron duality. The actual Th q is definitely non-analytic at any 
threshold (whether or not the amplitude vanishes at the threshold). Thus, the 
'radius of convergence' cannot correspond to the mass smaller than the threshold 
mass. Since in the actual QCD the thresholds exist at arbitrary high energy, the 
power expansion in Eq. (|186[ ) can be only asymptotic, with formally zero radius of 
convergence in I/tuq. 

In practice, the true threshold singularities are expected to be strongly sup- 
pressed at large energies, and the corresponding uncertainties in the OPE series 
quite small. Eventually they are expected to be exponentially suppressed, though, 
possibly, starting at larger energies. In the intermediate domain they can decrease 
as a certain power and must oscillate. This reflects the peculiarity of the Minkowski 
world. The terms left aside by the 'practical' OPE are exponentially suppressed in 
the Euclidean domain, ~ exp where Q 2 generically denotes the square of the 
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large momentum scale. Continuing this expression to Minkowski domain and taking 
discontinuity to determine transition probabilities, as discussed in Sects. 5 and 6, 
we observe oscillating rather than exponentially suppressed effects. 

As was illustrated in Ref . J74J , the power expansions like Eq. (|186|) are meaningful 



even beyond the power suppression where the duality-violating oscillations show up. 
In the case of the heavy quark widths where mass mq cannot be varied in experi- 
ment, the size of the duality- violating component may set the practical bound for 
calculating the widths. Thus, it is important to have an idea about its size. At the 
same time one should always include the leading QCD effects to the partonic expres- 
sions, rather than compare the actual observable with the bare quark result. In the 
model considered in Ref. |74[ , incorporating the power corrections from the practical 
OPE suppressed the apparent deviations by more than an order of magnitude. 

It is useful to keep in mind that violations of local duality, although so far poorly 
understood dynamically, are not arbitrary and must obey constraints following from 
the OPE. They cannot be blamed, for instance, for the systematic excess or system- 
atic depletion of decay probabilities; the actual width can only oscillate around the 
OPE predictions as a function of energy or quark mass, or the difference has to fade 
out exponentially. 

Although the simplest illustration of the asymptotic nature of the decay width 
1/rriQ expansion and related violations of local duality given above follows from 
the presence of hadronic thresholds, violation of local duality is a more universal 
phenomenon that is not directly related to existence of hadronic resonances nor 
even confinement itself. This has been illustrated in Ref. []67| by the example of 
soft instanton effects that do not lead, at least at small instanton density, to quark 
confinement - but do indeed generate computable oscillating duality-violating con- 
tributions to the total decay rates. 

One can parallel the OPE with expanding the interaction between quarks and 
gluons (and their related propagation) at small distances set up by the inverse energy 
scale in the problem. This is transparent in toy quantum mechanical problems 
where the terms in the OPE series can be traced back to the Taylor expansion of 
the potential V(x) at x^O. The behavior of V(x) in the vicinity of x = does not 
exhaust the problem, however. The spectrum of final states crucially depends on the 
asymptotics of the potential at large x; this is analogous to the resonances affecting 
local duality. There are more intriguing mechanisms associated with the finite- 
distance singularities of QCD interactions - the example is provided by instantons. 
The interested reader can find the inspiring discussion in the dedicated paper by 
Shifman in this volume [ |113 |. 



While conceptual grounds for violation of local duality have become more clear, 
its dynamical origin in actual QCD is still not well understood. What type of physics 
lies behind the finite- X singularities if they are relevant? How to quantify effects of 
resonances at high energies? We do not have answers yet. The instanton-based 
model |J7|, for example, while capturing correctly the gross features, clearly falls 
short in describing the size of the effects, at least under the standard assumptions. 
In all considered cases, actually, the effects of local duality violation turn out very 
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small in the asymptotic domain of large energies or quark masses, much below 
phenomenologically significant level. For practical purposes, we would like to know 
the limitations imposed by duality in the domain of intermediate energies where it 
a priori can be sizeable. The use of analytic methods can be limited here. 

Trying to get insights into the possible magnitude of violations of local duality 
in semileptonic decays of heavy quarks, recent paper 0] studied their numerical 
significance in the exactly solvable two-dimensional 't Hooft model. The model has 
built-in "hard" confinement related to linear Coulomb potential in 1+1 dimensions. 
Its spectrum consists of towers of infinitely narrow resonances. As was mentioned 
above, duality violation is not necessary related to existence of resonances. Never- 
theless, the intuition remains that resonance dominance is not "favorable" for the 
OPE, and problems might show up, for instance, through a delayed numerical on- 
set of duality, in that the approximate equality of the OPE predictions and the 
actual decay widths may set in only after a significant number of thresholds has 
been passed. To address such issues, the 't Hooft model seems to represent the most 
certain testing ground for local duality in the domain of decays of moderately heavy 
quarks. 

Contrary to naive expectations, surprisingly accurate duality was found between 
the (truncated) OPE series for T s \ and the actual decay widths. The deviations were 
suppressed to a very high degree almost immediately after the threshold for the first 
excited final state hadron is passed. No suspected delay in the onset of duality was 
found. Remarkably, the 'practical' OPE turns out to be the most efficient way to 
obtain numerical predictions for the total semileptonic width in the model to a high 
accuracy hardly accessible to direct computations, down to rather low quark masses. 

Some of the duality-violating features observed in those studies have natural 
explanations. At fixed energy release m^ — m c the magnitude of the deviations is 
smaller if mj, m c are both large than if they are both small. This is expected, since 
in the former case the heavy quark symmetry for the elastic amplitude additionally 
enforces approximate duality even when no expansion in large energy release can 
be applied. However, at fixed m& the duality violation decreased rapidly as m c 
decreases, in full accord with the OPE where the higher order terms are generally 
suppressed by powers of l/(mb — m c ). This is clearly a dynamical feature that 
goes beyond heavy quark symmetry per se, the quality of which deteriorates as m c 
decreases. 

To the extent the numerical findings of Ref. |5(| can be transferred to real QCD, 
violation of local duality in the total semileptonic widths of B mesons should not 
be an issue. The scale of duality violation lied far below the phenomenologically 
accessible limits, and could not affect the credibility of extracting |V^| or \V u b\. 

In reality there are, of course, essential differences between the two theories, 
including those aspects that are expected to be important for local duality (for 
a discussion, see Ref. j66|). Although many seem to optimistically suggest that 
duality violation is more pronounced in the 't Hooft model than for actual heavy 
flavor hadrons, some differences may still work in the opposite direction. In D = 2 
there are no dynamical gluons, nor a chromomagnetic field that in D = 4 provides a 
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significant scale of nonperturbative effects in heavy flavor hadrons. Likewise, there 
is no spin in D = 2, and no corresponding P-wave excitations of the spin-| light 
degrees of freedom which play an important role in D = 4. 

Two-dimensional QCD neither has long perturbative "tails" of actual strong 
interactions suppressed weakly (by only powers of logs of the energy scale). In 
D = 2 the perturbative corrections are generally power-suppressed, as follows from 
the dimension of the gauge coupling. As discussed in Ref. |66 |, it is conceivable 



that the characteristic mass scale for freezing out the transverse gluonic degrees of 
freedom is higher than in the "valence" quark channels. This would imply a possibly 
higher scale for onset of duality in a s /n corrections to various observables. 

Regardless of these differences, it was demonstrated that presence of resonance 
structure per se is not an obstacle for fine local quark-hadron duality tested in the 
context of the OPE. In the 't Hooft model resonances themselves do not demand a 
larger duality interval. As soon as the mass scale of the states saturating the sum 
rules in a particular channel (quark or hybrid) has been passed, the decay width can 
be well approximated numerically by the expansion stemming from the OPE. 

The same pattern is expected in QCD. To get a better idea regarding the onset of 
local duality, at least in semileptonic decays, it is important to study the saturation 
of the heavy quark sum rules, in particular the static SV sum rules discussed in 
Sect. 4. Are the values of /i^ and //g we use reached at the energy scale of 0.7 to 
1 GeV and right above we can use the perturbative description to account for higher 
states, - or we need to go higher in energy for that? Experimental study of actual 
b — > c transitions can give the answer. There are ideas how this can be studied on 
the lattices. The question of local duality in nonleptonic decays of beauty particles 
at present remains largely unknown. 



9 Conclusions and Outlook 

Heavy quark theory is probably one of the youngest branches of the QCD tree, yet 
it has become mature. Many practically important problems that laid dormant for 
years are now tractable. The development of the heavy quark theory, often initiated 
in the quest for the most accurate extraction of electroweak properties of quarks, 
enriched the library of available methods in QCD itself. At the same time, many 
natural limitations of QCD show up in the heavy quark expansion as well. Quark 
confinement is not completely understood yet, the infrared part of dynamics is more 
parameterized than solved. Therefore, every new result relying on the general prop- 
erties of the QCD interactions brings in a theoretical value. The success of the last 
decade was mainly related to incorporating Wilson's approach to the heavy quark 
expansion. It allowed to quantify the deviations from the heavy quark symmetry in 
many important cases, and led to a number of accurate predictions not appealing to 
the symmetry itself. The theory of preasymptotic corrections to the decay widths is 
one of such assets, whether or not b quark is heavy enough for precision applications 
of the heavy quark expansion. 
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A clear-cut recent manifestation of the power of the heavy quark theory is the 
framework it provides for determination of \V c b\. It is remarkable that the values 
of \V c b\ that emerged from exploiting two theoretically complementary approaches 
relying on quite different assumptions, are very close. The progress was not for 
free: it became possible only due to essential refinements of the available theoretical 
tools. They prompted us, in particular, that the zero recoil B^D* formfactor Fp* 
is probably closer to 0.9, with much larger deviations from the symmetry limit than 
previously expected. The predicted decrease in F D * brought more accurate experi- 
mental data which appeared later in a good agreement with \V&\ determined from 
T S \(B). The confidence in the latter method, in turn required many improvements 
beyond theory of nonperturbative effects as well, including accurate evaluation of the 
perturbative corrections. The latter story, although eventually successful, was far 
from simple (for the historical perspective see, e.g. Ref. [[J and references therein). 
Sharpening the numerical predictions was not possible without clarifying the nature 
of the heavy quark mass entering theoretical expression, the questions which long 
caused confusion in the literature. It is encouraging that now there is a consensus 
regarding also the numerical value of the running b quark mass - it is extracted by 
different groups from the T sum rules using the state of the art NNLO resummation 
of the perturbative effects. At the same time, a consistent theoretical definition of 
such a low-scale mass required investigating the dipole radiation in the non-Abelian 
gauge theory, a very general - although still somewhat theoretical - physical phe- 
nomenon. 

There are still problems to solve, ahead. Some of them are notoriously difficult, 
like violations of local duality. Others may look at once only technical, for example, 
developing the constructive Wilsonian OPE including perturbation theory, where 
the soft parts have to be removed from all diagrams preserving such basic properties 
of QCD interactions as gauge invariance. Like in the previous years, we may expect 
new links in the chain of theoretical advances to emerge from both types of studies. 
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